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Titla:  Conditional  Second  Order  Closure  For  Turbulent  Shear  Flows. 

Principal 

Investigator:  W.  Kollmann,  UC  Davis 

Summary 

^  The  research  work  during  the  first  year  was  concentrated  on  two 
areas:  The  foundation  of  conditional  closure  schemes  in  terms  of 

probability  density  functions  (pdf)  and  the  development  of  a  second  order 
closure  including  intermittency  factor  and  conditional  moments.  In  the 
theoretical  part  dynamics  of  single  and  multi-point  pdf's  for  velocity  and 
a  scalar  variable,  that  can  be  used  for  distinction  between  turbulent  and 
nonturbulent  zones,  were  considered  and  methods  of  closure  were 
investigated.  The  transport  of  apparent  stress  in  the  nonturbulent  zone 
of  turbulent  shear  flows  with  a  free  boundary  was  included  in  the  second 
order  model  in  terms  of  their  dynamic  equations.  Conditions  governing  the 
effect  of  the  fluctuating  interface  on  mean  velocity  and  apparent  stress 
in  both  zones  were  established  and  closure  models  were  put  forward.  The 
resulting  second  order  closure  was  compared  with  experiments  for  several 
plane  shear  flows  and  good  agreement  was  found. 

Research  Objectives 

The  objective  of  the  proposed  research  project  is  the  development  of  a 
second  order  closure  model  for  conditional  moments  and  the  intermittency 
factor.  The  foundation  of  the  closure  scheme  are  to  be  investigated  and 
the  resulting  model  should  be  applicable  to  a  wide  range  of  turbulent 
shear  flows  with  free  boundaries. 


The  research  work  on  this  project  started  in  July  1984  with  M 
Hortazavi  as  a  Ph.D.  student  and  S.  Byggstoyl  fro  TU  Trondheist  (Borway)  as 
a  Postdoctoral  Fellow. 

M.  Hortazavi  and  W.  Kollmann  worked  during  the  first  year  on  the 
probability  density  formulation  (pdf)  of  conditional  closures.  M. 
Hortazavi  started  with  a  thorough  survey  of  the  existing  literature  on 
pdf-methods.  This  lead  to  the  definition  of  his  contribution  as  the 
single  and  multi-point  pdf  formulation  of  conditional  closure  schemes.  He 
then  began  working  on  the  specific  problem  of  turbulent  diffusion  of 
prassive  scalars,  for  which  he  is  currently  developing  closure  models 
based  on  the  work  of  Lundgren  (1],  (2],  Ievlev  [3],  Kuo  and  O'Brien  [4], 
Kollmann  and  Janicka  [5]. 

S.  Byggstoyl  and  W.  Kollmann  developed  the  second  order  closure  [6] 
further  by  considering  the  transport  of  apparent  stress  in  the 
nonturbulent  zones  of  shear  flows  with  free  boundaries.  It  was  found  that 
conditional  mean  velocities  and  conditional  stresses  undergo  additional 
transport  and  production/destruction  due  to  the  random  fluctuations  of  the 
interface  separating  the  turbulent  and  nonturbulent  zones.  This  effect  is 
represented  by  additional  terms  in  the  transport  equations  for  these 
moments  and  it  was  shown  that  those  terms  were  linked  by  a  consistency 
condition  and  approach  particular  limit  forms  as  the  distance  from  the 
turbulent  shear  layer  increases  (Corrsin-  Kistler  relations).  These 
theoretical  results  are  contained  in  appendix  I.  The  second  order  closure 
model  was  then  complemented  with  the  transport  equations  for  the  apparent 
stress  in  the  nonturbulent  zone.  The  condition  that  the  fluctuations  in 


th*  non turbulent  son*  *r*  irrotational  loads  to  decay  lairs  with  distanca 
frosi  th*  turbulant  ration,  which,  together  with  tha  consistancy  condition, 
rastriet  th*  possibl*  closura  expressions.  Th*  numerical  solution  of  tha 
resulting  system  of  nonlinear  parabolic  differential  aquations  required 
simultaneous  solution  of  th*  stress  aquations  in  each  son*  in  order  to 
cop*  with  nonlinear  instabilities  which  frequently  occurred  in  th* 
sequential  solution  SMthod.  Thus  a  blocksolver  was  introduced  for  the 
stress  equations  which  improved  th*  stability  characteristics 
signif  iciently  and  lead  to  a  sod*  rate  gain  in  computing  time .  Th* 
comparison  of  th*  results  with  experiments  in  several  plan*  shear  layers 
is  quite  satisfactory.  The  results  are  presented  in  detail  in  appendix  II. 

Th*  research  on  conditional  closure  schesss  lead  to  a  new  idea. 
Multi-scale  models  based  on  conditional  statistics  using  multizonal 
distinction  with  respect  to  an  appropriate  scalar  variable.  This  aspect 
of  the  project  is  discussed  in  detail  in  appendix  III  which  contains  also 
th*  same  preliminary  results. 
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Appendix  I:  Stress  transport  in  the  rotational  and  irrotational 

zones  of  turbulent  shear  flows 

S.  Byggstoyl,  W.  Kollmann 
Department  of  Mechanical  Engineering 
University  of  California 
Davis,  California  95616 

AIP  Classification:  47.25.Fj,  47.25.Jn 

Abstract:  The  transport  equations  for  intermittency  factor  and 

conditioned  moments  are  analyzed  for  turbulent  shear  flows  with  free 
boundaries.  Conditions  are  established  for  molecular  diffusion  to 
dominate  the  progression  of  the  turbulent  zone.  The  terms  related 
to  the  dynamics  of  the  interface  for  mean  velocity  and  apparent 
stresses  in  turbulent  and  nonturbulent  zones  are  shown  to  be  linked 
by  a  local  relation  and  the  limit  of  large  distance  from  the 
turbulent  region  for  the  nonturbulent  zone  stress  equation  is  given. 
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Stress  transport  In  the  rotational  and  irrotational 


zones  of  turbulent  shear  flows 

S.  Byggstoyl,  W.  Kollmann 
Department  of  Mechanical  Engineering 
University  of  California 
Davis,  California  95616 


1.  Introduction 

Turbulent  shear  flows  with  free  boundaries  such  as  jets  and  mixing 
layers  show  rotational  and  irrotational  fluctuations  of  velocity  near  the 
free  boundary.  This  was  established  by  Corrsin  and  Kistler  [1],  [2] 
experimentally  for  several  boundary-layer-type  flows.  Since  then  a  large 
number  of  experimental  results  became  available  (see  references  given  in 
[3]).  The  theoretical  investigation  of  intermittently  turbulent  flows 
started  with  Corrsin  and  Kistler  12],  Phillips  [4],  17]  Stewart  [5], 
Corrsin  and  Phillips  (6],  dealing  in  particular  with  interface  properties 
and  irrotational  fluctuations.  Lumley  [8],  [9]  introduced  the  statistics 
of  multi-valued  random  functions  from  the  treatment  of  interface 

fluctuations  and  developed  closure  ideas.  The  first  closure  model  based 

on  first  and  second  order  moments  was  published  by  Libby  [10],  [11]. 

Dopazo  [12]  and  Dopazo  and  O'Brien  [13]  established  the  theoretical 

representation  of  interface  related  processes  in  conditioned  moment 
equations.  Pope  [14]  showed,  that  conditioning  can  be  easily  applied  to 

pdf -equations . 

The  present  paper  extends  the  analysis  of  ref.  [3]  to  stress  transport 
in  the  nonturbulent  zone  of  shear  flows  with  free  boundaries.  The 


condition  of  irrotationality  in  terms  of  the  Corrsin-Kistler  equation  is 
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The  statistical  description  of  turbulent  flows  can  be  refined  by 
conditioning  of  expectations  to  capture  particular  properties  of  the  flow 
in  an  explicit  fashion  [3], [12].  For  this  purpose  a  non-negative  scalar 
variable  is  selected  in  such  a  way,  that  a  local  condition 
being  satisfied  at  a  point  ^£,4)  of  the  flow  field  corresponds  to 


>  h  >  O 

where  Ht  is  the  threshold  value,  and 


<h 

corresponds  to  the  condition  being  violated.  In  the  present  case  the 
condition  is  that  the  flow  is  turbulent  at  the  point  considered  and  thus 


>0 


is  taken  as  discriminating  variable.  It  should  be  noted  that  many  other 
conditions  can  be  considered  (such  as  hot-cold,  bumt-unbumt, 
colored-clear)  which  would  lead  to  different  variables  ) 


Conditioning  of  flow  variables  can  be  done  with  the  aid  of  the  indicator 
function 


for  4)  >  h 


ottier*stse 


The  locations  X  where 


S/'X)  i  )  m  -  h  m  O 


(!) 


form  a  surface  separating  turbulent  from  non- turbulent  zones.  This 

§ 

surface  propagates  with  velocity  Z ^  and  progresses  relative  to  the  fluid 
in  its  normal  direction  (positive  into  turbulent  zone)  with 

speed  1/  p2] 

-  %  -  Vs  (2) 


Then  follows  (3) 

§f  -  - 

and 

2 I  V*  , 

Where 


n) 


to) 


/F5/f?5) 


Conditional  momenta  can  now  be  defined  for  the  turbulent  cone 


y>*  * 


and  the  non-turbulent  zone 

*  '~t 


9 -? 


where  the  intermitt ancy  factor  is  given  by 


rr) 


r  ■  <T> 


Points  on  the  interface  move  with  velocity  ZJ,5  thus 


SI  + 

96  *  V«  ^ 


O 


which  can  be  recast  in  terms  of  fluid  velocity 
progression  velocity  of  the  interface. 

91 


^  +  V«KL 


9z 

£ 


(*) 


H) 


and  the  relative 


no) 


Averaging  leads  to  the  exact  equation  for  the  intermittency  factor  (3} 


k  *  “  < 'rJh>> 

Introducing  the  unconditional  mean  as  convection  velocity  yields 

k  +  4  =  'A  ~ 71 

This  form  of  the  transport  equation  for  the  intermittency  factor  shows. 


that  the  turbulent  diffusion  of  y~  i*  due  to  the  relative  movement  of 
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turbulent  and  non-turbulant  zones  if  a  point  moving  with  the  unconditional 
mean  velocity  is  followed.  The  intermittency  source  S ^ 


ns) 


is  the  rete  at  which  the  volume  of  the  zone  (f>  >  h.  grows  per  unit  volume 
of  fluid.  A  detailed  analysis  (3}  shows  that  S^.  can  be  represented  in 
several  fonts.  The  dependence  of  this  growth  rate  on  the  scalar 
variable  discriminating  between  the  zones  can  be  shown  explicitly  (3). 
If  <f>  satisfies 


then  f o 1 lows  for 


m) 


where  2>/Pi  denotes  the  substantial  derivative.  The  transport  equation 
for  (p  yields  then  the  representation 


sr  -  n4n§£)  +  o]xs>> 

This  result  shows  that  two  mechanisms  for  growth  are  present  namely 
diffusion  and  sources.  The  influence  of  the  source  Q  on  depends 

crucially  on  the  threshold  level  H  and  the  limit  of  Q  as  (fi 
approaches  zero.  If 

CGJ(S)>  =  C(h*) 

h  -* o 

with  an  exponent  /  ,  then  is  the  effect  of  { C'^TS) ^ 

on  Sjft  negligible  if  H  is  sufficiently  small.  For  this  case  can 


VW V  v‘V.V V  ^ V*  s'.  ^ 


A -  .  A  .‘i^YA.V.t  '.tWJ 


be  concluded,  that  discrimination  with  a  conserved  scalar  yields  the  same 
result  as  discrimination  with  a  non-conserved  scalar.  This  does  not  hold, 
if  the  threshold  H.  is  raised  to  a  level  with  large  rates  of  production 
or  destruction  of  <$>  An  estimate  for  the  threshold  level  h  ,  such 
that  the  diffusive  part  of  S^.  dominates,  can  be  obtained,  if  the  source 
term  Q  can  be  represented  as  power  in 

Set)  =  *  <P  n 


for  small  <p  The  source  term  S y  is 


r 

^  £<5? 


7*  t 

'3 


where  Sj  and  can  be  written  in  terms  of  conditional  expectations 


and 


S/  2  ( 'q/0 -h)  PfA) 


it  6?  is  a  power  of  tfi  we  obtain 


s? 


A  h  *  P(A) 


The  diffusive  part  can  be  estimated  by 


Ofsf)-  Pek)fji 


Ob) 


07) 


The  length  scale  /  can  be  related  to  the  macro-scale  L  of  the  flow,  if 


it  is  required,  that  the  dominant  production  and  destruction  terms  in 
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the  equation  for  scalar  dissipation  have  the  same  order  of  magnitude. 


Then  follows 


oft)  -&)'* 


±  sr 


and  the  relative  order  of  magnitude  of  o^,  to  is 


o(-$) 


Where  Z.  and  U  are  the  macro-scales  for  length  and  velocity 


and  7^  —  uL/\>  .  If  this  ratio  is  required  to  be  less  than  unity,  we  get 


t.  tl“/  *  c  t/22  j  t-LL 
h  <  Q  (tj  LA 


Thus  the  following  conclusions  are  reached.  If  £}f<f>)  approaches  zero 
faster  than  (f>  ,  the  diffusive  part  of  Sp  is  always  dominant  for 


sufficiently  large  Re-number.  If  St  is  less  than  unity  however,  we  find 


*-*>**&)*& 


and  the  threshold  level  H.  for  Cp  is  restricted  by  the  Re-number  from 
below.  If  hi  violates  this  inequality  the  source  term  will 


dominate  the  development  of  the  intermittency  factor. 


The  structure  of  the  diffusive  contribution  of  Sy.  can  be  illustrated 


for  the  case  of  Gaussian  statistics  of  the  scalar  (f>  This  example  can 
be  viewed  as  Cf)  being  lognormal  or  i* 1 4>  being  an  unbounded 
discriminator.  The  mean  and  the  variance  are  taken  as 


constants  and  the  spatial  correlation  coefficient 


+  ho  i- 
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is  developed  in  a  Taylor  series.  Writing  for  the  source  term 

)J«p  -  h j)  =  £^o 

■+  -  *  frfSr*  -*»/ 


|  and  introducing  the  assumption  that  and  )  are 

» 

■  Caussian-distributed ,  we  get  after  some  algebra 


I 


~h>r£^/0)  m 

This  result  has  several  interesting  properties.  First  it  shows  that  the 
intermittency  source  depends  on  the  probability  for  (fif* )  ~  ft-  at  the 
point  X  and  therefore  on  the  threshold  value  /?.  for  the  discriminating 
scalar  Furthermore  it  becomes  evident  from  (19),  that  for  the  limit 
of  infinite  Re/Sc-numbers  the  intermittency  source  will  approach  a  nonzero 
and  bounded  constant,  because 


> 

I 

! 

I 


ro) 


becomes  independent  of  Re/Sc-numbers  by  virtue  of  the  same  arguments  that 
apply  to  viscous  dissipation.  Hence  can  (19)  be  recast  as  follows  for  the 
example 


Sf  -  fflf 


showing  explicitly  the  dependence  on  the  scalar  time  scale.  The  influence 
on  the  velocity  fluctuations  enters  via  the  pdf  ^  of  the  scalar  « 


9 


which  is  transported  and  diffused  by  the  velocity  field.  Finally,  it 
should  be  noted  that  the  intend ttency  source  is  not  positive  definite, 
but  can  become  negative  for  f $>)  <  h. 


3.  Mean  velocities  in  turbulent  and  non-turbulent  rones. 

The  balance  of  momentum  can  be  averaged  conditionally  thus  leading  to 

the  equations  for  tone-averaged  mean  velocities.  The  turbulent  zone  mean 

jct 

velocity  satisfies 


and  for  the  non-turbulent  zone  follows 


,7  54 


az) 


Both  equations  contain  a  term  representing  the  effect  of  the  interface 
movement  on  the  average  momentum  in  the  two  zones.  These  interface  terms 
can  be  given  in  the  form 


4'  * 

and 

ST'  •  '%nAtfs>> 


(23) 


(2*) 


where  /  m  *  for  the  turbulent  con* 
non- turbulent  zone.  The  interface  tenu 
independent  but  linked  due  to 


quantities  and  /  «  O  for  the 

^  0  ^  (5  f  1  ft  y  O 

^  \  mnA  «•  not 


%  -/4  * 

Combination  of  (21)  with  (22)  must  reproduce  the  equation  for  the 
unconditioned  mean  velocity.  Thua  follows  [16]  after  game  manipulations 


DL 


/4  -  <2  )f'sf  -  **> 


and 

f?S)  >  -  ^  at>) 

Iquations  (25)  and  (26)  must  be  satisfied  for  closure  expressions 
relating  7^'  and  2$))  to  other  quantities  and  are  called 

consistency  conditions. 

So  far  no  use  has  been  made  of  the  fact  that  the  non- turbulent  zone  is 
irrotational .  Following  Corrsin  and  Kistler  [2]  it  is  required 


At £/«  Ofh*)  for  -T  ~  O 


or 


(27) 


where  hi  is  the  threshold  value  for  the  turbu  lent -nonturbu  lent 


discrimination.  This  relation  can  be  applied  to  establish  the  properties 
of  transport  equations  for  moments  in  the  non- turbulent  zone  as  the  outer 


11 


edge  of  the  flow  is  approached  (corresponding  to  approaching  zero). 

Averaging  (27)  laads  to  the  Corrsin-Kistler  equation 

°  'fit'll  ~  &t)  *  ^ -  VH**  )Jh)>  -  0  #*) 

of  first  order.  The  interface  term  in  (28)  can  readily  be  shown  to 
satisfy  the  following  conditions: 

-o 

o 

and 

- 0 

The  first  limit  follows  from  the  fact  that  in  the  non-turbulent  flow  field 
the  mean  vorticity  approaches  zero  as  the  distance  from  the  turbulent  zone 
grows.  The  second  limit  follows  directly  from  (28).  Therefore  is  a 
representation  of  the  interface  term  in  (28)  in  the  form 

~  ^'•'1  V&>  -  a,) 

possible  with  an  unknown  bounded  function  These  properties  of  the 

■$o 

Corrsin-Kistler  equation  will  be  used  in  the  discussion  of  the  dynamics  of 
the  apparent  stress  tensors  in  turbulent  and  non-turbulent  zones. 

Corrsin-Kistler  equations  of  higher  order  can  be  obtained  from  (27)  by 
multiplication  with  fluctuating  components  and  averaging.  Thus  follows 

tec) 


Where  A  °  2  JJ  e  fJ  ° 

2  <x  n 


denotes  the  kinetic  energy  in  the 


'  ■  v'.' V.V/.' .v.v v.  v  ■.  .v. 


.•  /.  .v  .• 
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non- turbulent  zone.  For  the  limit  j/  — ?  O 

3 


the  relation  [2], [5] 


^  -2r 

_  £>k° 


**  *'**  * 


is  obtained,  which  shows  that  the  effect  of  the  apparent  stress  on  the 
mean  velocity  in  the  non-turbulent  zone  becomes  analogous  to  the 
pressure-gradient . 


4.  Apparent  stresses  in  turbulent  and  non-turbulent  zones. 

The  momentum  balances  and  mass  conservation  can  be  multiplied  with 
appropriate  fluctuations  and  the  indicator  function  and  averaged  to 
establish  the  transport  equations  for  the  zonal  stress  tensors.  Thus  the 
equation  for  the  stress  tensor  in  the  turbulent  zone 
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a/) 


is  obtained  in  the  form 


2J* v* 
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-is* 
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The  tensor  of  apparent  stresses  in  the  non-turbulent  zone 
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S3S) 


is  governed  by 
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The  substantial  derivative  in  (32)  and  (34)  is  defined  by 

-  -£  y.  A)  ^2. 

3>ir  ~  $6  *  ^3*. 

where  the  superscript  applies  to  the  velocity. 

Both  stress  equations  (32)  and  (34)  contain  dissipation  rates 
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and 
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Their  importance  depends  on  the  discriminating  scalar.  For  the  present 
case  of  distinction  between  turbulent  and  non-turbulent  zones  becomes  the 
dissipation  in  the  non-turbulent  zone  negligible  and  the 
unconditioned  rate  of  dissipation  is  related  to  the  conditioned  rates  by 


If  the  discriminating  scalar  is  chosen  to  afford  distinction  between  other 


properties  of  the  flow  such  as  hot  and  cold,  then  both  zones  can  be 
turbulent  and  both  dissipation  rates  and  are  important. 


The  turbulent  fluxes  in  (32)  and  (34  are  defined  by 


S*  = 


-  vyrv/ 


-  °L^p' 


-  eL 
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and 
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The  existence  of  a  fluctuating  interface  separating  turbulent  and 
non-turbulent  zones  leads  to  the  transport  of  apparent  stress  through  the 
interface,  transport  due  to  interface  movement,  production  of  apparent 
stress  due  to  interface  fluctuations.  These  effects  are  represented  by 
the  term  group  defined  by 


%m 


f  f/3  )J?S)  >  +  P)  ,  ° 
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and 
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The  interface  terms  S*  and  S°  are  not  independent.  The  relation 

e<o 

between  conditioned  and  unconditioned  stress  tensors 

£7 Jig  -  f  {gTj*  +  //-  f) ~Vl} 


the 


leads  to  a  consistency  condition  for  the  Introducing 

abbreviations 

»  4  -  £ 

and 

23/)  ~p  -ft 

we  obtain  after  some  manipulations  using  (32),  (34),  (41)  the  following 
relation 

5^  -  sj  -  -g- 

-  //  444  $')  y-  /v  -/  ;  z^,  ) 

•*  +AV*§i) 

£  V  V  ¥  ¥ 
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This  rather  complicated  relation  can  be  somewhat  simplified  if  the  limit 
case  of  high  Re-numbers  is  considered.  Then  follows  that  the  viscous 
terms  in  (42)  can  be  neglected,  because  they  scale  with  the  mean  fields. 


The  transport  of  apparent  stress  in  the  zone,  where  the  discriminating 
scalar  0  is  below  the  threshold  level,  merits  discussion  in  particular 
if  distinction  between  turbulent  and  non-turbulent  zones  is  considered. 
We  consider  the  terms  on  the  right  hand  side  of  (34)  for  this  case  in 
turn.  The  production  of  stress  due  to  the  interaction  of  stresses  and 
mean  deformation  rates 


—  7J7J 


can  be  analyzed  using  the  Corrsin-Kistler  relation  (28).  Taking  (29)  into 
account  we  find 

jSf  — ’  A  P*  /v 


£  “  rffc  ~  V*V'°  ix?  ^  V  Ku  +  VWKr  ) 

Hence  we  find  for  the  limit 
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On  the  other  hand  is  for  the  limit  f4'  — *  /  the  original  form 
of  &  relevant  because  the  non-turbulent  patches  become  increasingly 
rare  in  this  limit  and  will  follow  the  motion  of  the  surrounding  turbulent 
fluid.  It  is  instructive  to  write  out  the  production  terms  for  both 
limits  in  case  of  a  plane  parabolic  flow  (i.e.,  plane  jet). 
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Mote  that  for  parabolic  flows  derivatives  with  respect  to  X')  ,  can  be 

neglected  (boundary  layer  assumption).  It  becomes  apparent  from  the  table 
above,  that  energy  is  fed  into  the  stress  component  in  the 

limit  jj*  — *  O  ,  whereas  for  ^  /  2J,  1  receives  the  energy  as 

in  the  turbulent  zone.  The  limit  y-  — *  Q  is  consistent  with  the 


relation  among  the  normal  stresses  (valid  for 
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— >  O  ) 


**  =  V, ol 


vcz 


(Ul) 


in  the  nonturbulent  zone  which  was  obtained  by  Phillip  [4]  and  Stewart  [5]. 

The  diffusive  flux  contains  the  triple  correlations  of 
velocity,  which  can  be  analyzed  by  means  of  the  Corrsin-Kistler  equation 
for  triple  moments 


+  2(/-r)  --  -  z±°2<\-  vA°Jc\)ct*rs>y 

+  OCh**-)  ft*) 

where  ■/ J  denotes  the  rate  of  strain. 

For  the  limit  ^  >  Q  the  relation 


is  obtained.  Hence  is  the  divergence  of  the  flux  of  apparent  stress  for 
this  limit  composed  of  gradients  of  the  flux  of  kinetic  energy  and  the 


correlation  of  kinetic  energy  with  the  fluctuating  strain  rate.  More 
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distance  from  the  interface  increases.  Phillips  (71  showed  that  for 
large  y 

OfvJ)  =  y~2 

Where  y  is  the  normal  distance  from  the  interface.  Using  the  Bernoulli 
equation,  which  is  valid  in  the  non-turbulent  zone  only,  it  follows  that 

0{p°)  *  y 


These  estimates  can  be  used  to  assess  the  relative  importance  of  several 
terms  in  the  moment  equations  for  the  non-turbulent  zone  for  large 
distance  from  the  turbulent  zone  or  as  y  —>  O  It  is  worth  noting 
that  irrotational  fluctuations  are  induced  by  boundaries.  Thus  is 
statistical  homogeneity  not  possible  and  decay  estimates  have  to  be  used 
instead. 

The  interface  terms  in  (39)  show  that  the  growth  of  the  turbulent  zone 
leads  to  an  increase  of  non-turbulent  zone  stresses  by 


/ 


,-f  S,W 

Which  may  be  cancelled  by  the  complex  production/destruction  group  of 


terms 


,  as  the  consistency  relation  (42)  indicates,  which  contains 


the  intermittency  and  conditioned  momentum  sources. 


5.  Conclusions. 

The  transport  equations  for  intermittency  factor  and  conditioned 
moments  were  set  up  and  their  properties  were  analyzed.  The  conclusions 


can  be  summarized  as  follows. 

(1)  Intermittency  factor:  The  intermittency  source  was  shown  to  be 

coaposed  of  two  different  terms  representing  growth  due  to  molecular 
diffusion  and  production  of  the  scalar  variable  used  for 

discrimination  between  the  zones.  The  diffusive  term  is  the  dominant 
mechanism  for  growth  if  the  source  term  approaches  zero  faster  than 
the  first  power  of  the  scalar.  further  was  shown  that  the 
intermittency  source  is  not  always  positive. 

(2)  Conditioned  mean  velocity:  The  equations  for  turbulent  and 

non- turbulent  zone  mean  velocity  contain  source  terms  describing 
siomentum  transport  through  the  interface  and  production/destruction  of 
mean  momentum  due  to  interface  movement.  These  source  tins  are  not 
independent  but  linked  by  a  local  consistency  relation.  Which  shows 
that  the  difference  between  these  sources  is  proportional  to 
intermittency  source  and  the  difference  of  conditioned  mean 

velocities.  The  Corrsin-Kistler  equation  introduces  the  condition  of 
irrotationality  in  the  non- turbulent  zone  and  proves  that  for  the 
limit  £  O  the  effect  of  apparent  stress  on  the  non- turbulent 
mean  velocity  becomes  analogous  to  the  mean  pressure-gradient. 

(3)  Conditioned  stress  tensors:  The  transport  equations  for  the  apparent 
stress  tensors  in  turbulent  and  non-turbulent  zones  contain  source 
terms  describing  the  production/destruction  of  stress  due  to  the 
interface  fluctuations.  These  source  terms  are  not  independent  but 
linked  by  a  local  consistency  relation  as  in  case  of  the  mean 
velocities.  The  analysis  of  the  stress  transport  equation  in  the 
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non-turbulent  rone  shows  that  production  in  the  limit  — >  C>  for 

boundary- layer- type  conditions  shifts  from  to  ,  thus 

confirming  the  dominance  of  component  *  for  this  limit.  The 

Corrsin-Kistler  relation  for  triple  correlations  shows,  that  under  the 
same  conditions  becomes  diffusion  of  dominant,  if  the 

correlation  of  kinetic  energy  and  strain  rate  is  weak. 
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Appendix  II:  A  closure  model  for  conditioned  stress  equations 

and  its  application  to  turbulent  shear  flows 

S.  Byggstoyl,  W.  Kollmann 
Department  of  Mechanical  Engineering 
University  of  California 
Davis,  California  95616 

AIP  Classification:  47.25.Fj,  47.25.Jn 

Abstract:  A  second  order  closure  model  based  on  intermittency  factor  and 
conditioned  moments  is  developed.  The  transport  equations  for  the 
nonturbulent  zone  stresses  are  included  in  the  model.  The  resulting  model 
is  then  compared  with  measurements  in  several  shear  flows  and  satisfactory 
agreement  between  calculation  and  experiment  is  obtained. 


k  cloture  model  for  conditioned  stress  equations 


and  its  application  to  turbulent  shear  flows 

S.  Byggstoyl,  W.  Kollmann 
Department  of  Mechanical  Engineering 
University  of  California 
Davis,  California  95616 

(1)  Introduction 

Turbulent  shear  flows  with  free  boundaries  such  as  jets,  wakes,  mixing 
layers  show  an  intermittent  character  in  the  fluctuations  of  velocity 
changing  from  rotational  to  irrotational  and  vice  versa  (1],  The 
prediction  of  this  type  of  shear  flows  can  be  based  on  unconditional 
moments  [2], [3]  on  conditional  moments  [4], [5], 16].  The  latter  case  is 
considered  here.  Which  allows  calculation  of  the  intermittency  factor  and 
statistical  moments  characterizing  the  fluctuations  in  the  individual 
zones.  The  closure  model  developed  in  this  paper  is  based  on  the  work 
presented  in  [7], [8] ,19).  The  new  part  is  the  inclusion  of  transport 
equations  for  the  apparent  stresses  in  the  nonturbulent  zone.  The 
properties  of  turbulent  and  nonturbulent  zone  moments  discussed  in  [9]  are 
used  in  the  development  of  closure  expressions.  The  resulting  model  is 
then  compared  with  measurements  in  several  shear  flows  and  satisfactory 
agreement  between  calculation  and  experiment  is  obtained. 

(2)  Intermittency  factor 

The  equation  for  the  intermittency  factor  [5], [8] 

sf  +  +  sr  n) 


contains  tha  rata  of  growth  of  tha  turbulant  zona  at  the  axpanse  of 
tha  nonturbulant  zona.  In  raf.  18}, [9]  savaral  representations  of  -S^, 


Sy  3  fi/StS)) 


r 2 ) 


ware  discussed.  Vote  that  )/  denotes  the  relative  progression 
locity  of  the  interface  and  res)  is  defined  by 


va 
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with  S  0  being  the  implicit  equation  for  the  interface  (see 
(5), (8]).  Since  no  representation  of  5^  in  terms  of  first  and  second 
order  moments  only  is  available,  a  closure  model  is  required  in  the 
context  of  second  order  closures.  Several  closures  have  been  suggested 
[4], [7] ,[10]  for  the  source  Sj,  We  follow  here  the  model  suggested  in 
ref.  [71. 


-  o. 
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with  •*»  J.S  t  “  OSS'  t  O.OS~ .  The  first  part  in  (4) 

represents  growth  of  ^  due  to  production  of  apparent  stress  in  the 

turbulent  zone.  The  second  term  on  the  right  hand  side  of  (4)  reflects 
the  transport  of  mass  and  momentum  due  to  spatial  inhomogeneity  thus 
increasing  the  intermittency  factor.  The  last  term  in  (4)  destructive. 


with  downstream  distance.  Subsequent 


that  in  fact  decn 

distortion  of  the  wake  flow  lead  to  production  of  turbulence  and  immediate 
increase  of  intermittency  factor  [11].  Thus  is  expected  to  depend  on 
the  mean  strain  rate  as  suggested  in  (4). 

The  difference  of  the  mean  velocities  in  the  turbulent  and 
nonturbulent  zones  appearing  in  the  intermittency  equation  (1)  is  in  fact 


the  turbulent  diffusion  of  .  This  difference  for  the  cross-flow 

component  is  for  parabolic  flows  estimated  by 

4  -  4  -  ct  -g  4*2?  j ^  cq-  <?'*>> 

because  of  numerical  inaccuracy  in  the  calculation  of  the  cross-flow 

— —  <v 

component  of  2^  ,  ££  in  boundary- layer- type  flows.  For  elliptic 

flows  are  all  momentum  balances  included  in  the  system  of  equations  and 
(5)  would  be  avoided. 


[3) Mean  velocity  in  turbulent  and  nonturbulent  zones 


The  exact  equation  for  the  turbulent  zone  mean  velocity  ZJ. 


m  y 
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and  for  the  nonturbulent  zone  2^ 
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contain  interface  terms  and  They  are  linked  by  the 

consistency  condition  [7] ,[8] 

4'  -  4°  - 

for  high  Re-number  are  the  viscous  terms  in  (8)  negligible  and  we  obtain 

4'  -  4°  -  &-  &  4  *  *> 

Hence  only  one  of  the  momentum  sources  or  requires  closure. 

It  is  helpful  for  the  construction  of  a  closure  to  consider  the 
nonturbulent  zone  equation  (7).  The  apparent  stress  term  in  (7)  can  be 


recast  as  follows 


U’V?  A 


The  first  part  has  the  proper  divergence  form  of  a  diffusive  term  whereas 
the  second  part  is  a  source  term.  Comparing  (10)  with  the  corresponding 
expression  in  the  turbulent  zone  (eq.  (16) 

M  *  4 

a  similar  source  term  but  with  the  opposite  sign  appears.  Considering  a 
shear  flow  such  as  a  jet,  it  becomes  clear  that  for  the  turbulent  zone  is 
the  second  part  is  a  genuine  source,  whereas  for  the  nonturbulent  zone  it 
is  a  sink,  which  leads  to  wave-like  variations  of  if  included  in 


the  equation  (7).  Since  there  is  no  experimental  or  theoretical  evidence 
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for  such  solutions,  this  sink  in  (7)  must  be  balanced  by  a  corresponding 
term  in  5*  .  Furthermore  is  momentum  transfer  between  turbulent  and 
non  turbulent  zones  only  possible  in  the  mean  if  2^  and  2^  are  not 
equal.  Hence  the  following  closure  model  for  J?,0  is  suggested 

4'  -  r//) 

with  Cl.  ^  J-O  .  Hote  that  S°  vanishes,  if  2^  —  for  a  region 

y  *v 

with  finite  volume,  which  implies  that  approaches  zero  also. 

The  consistency  condition  (9)  yields  then 

4'  -  -  5  of  +  srJ 

+ f  ft  f,z) 

Hote  that  for  the  longitudinal  component  (etm  / )  the  pressure  term  in  (12) 
becomes  negligible.  The  apparent  stresses  in  both  zones  are  included  in 
the  system  of  variables  and  therefore  conclude  (11)  and  (12)  the  closure 
of  the  mean  velocity  equation. 


(4)  Stresses  in  turbulent  and  nonturbulent  zones 

The  exact  transport  equations  for  the  apparent  stresses  in  turbulent 
and  nonturbulent  zones  follow  from  mass  and  momentum  balances  [9],  They 
can  be  given  for  high  Re-numbers  in  the  form 
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The  diffusive  fluxes  are  defined  by  [9] 
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The  interface  terms  *  and  £  °  are  defin 

satisfy  the  consistency  relation  (see  Eq.  (42)  in  [9]) 
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Where  422,  =  <^  -  2^  • 


^  47 


This  relation  (17)  is  valid  for  high  Re-numbers  and  Ap>  ~  J£  **  O 

It  is  reasonable  (at  least  for  thin  shear  layers)  to  assume  that  the  mean 

^  — 

pressures  in  the  two  zones  are  equal  because  yD  “ * ^  for  ^  ~>  0 
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(outer  edge)  and  JQ  f°r  ^  /  (center  of  shear  layer)  and yC> 

is  constant  across  a  thin  shear  layers  in  first  order  approximation. 


Hence  is 


a  good  approximation  to  first  order. 


The  closure  of  the  stress  equations  is  concerned  with  two  groups  of 
terms:  The  classical  terms,  which  are  counterparts  of  terms  in  the 
unconditional  stress  equations,  and  the  interface  terms.  For  the  closure 
of  the  classical  term  group  existing  models  13], [12], [7]  are  carried  over 


f  rom  the 


unconditional  to  the  conditional  correlations.  This  is 


certainly  justified,  because  closure  assumptions  based  on  quasi-Gauss ian 
behaviour  of  higher  moments  or  relaxation  to  Gaussianity  in  the  absence  of 
strain  rates  and  boundaries  (such  as  Lumley's  diffusion  closure  [12])  hold 
with  better  accuracy  for  conditioned  moments  than  for  unconditioned 
moments.  This  follows  from  the  observation,  that  conditioning  removes  the 
spike  in  the  pdf  corresponding  to  the  other  (nonturbulent)  zone  thus 
bringing  in  particular  flatness  factors  closer  to  the  Gaussian  values. 
This  is  backed  up  by  experimental  results  [13], [14].  The  interface  terms 
on  the  other  hand  require  new  considerations. 


4.1.  Closure  of  the  turbulent-zone-stress  equations 

The  closure  of  the  classical  terms  requires  only  brief  discussion, 
because  established  closure  models  will  be  modified  for  conditioned 
correlations.  The  dissipation  of  stress  is  taken  in  high  Re-number  form 


•^4,  3 


where  the  dissipation  rate  £?  is  given  by 


e  s  rf  e 

z  — 


8 


For  the  diffusive  flux 


correlations  is  applied 


Lumley's  model  [12]  for  the  triple 


where 


iu -  vtv;  +  V?vj£  v:v/  +  WW&t&j!  ao) 


“7 


4  *  *  ^*4  ^ 


The  constant  Cj_  was  set  to  <^  =  7. 5”  .  The  pressure  correlations 
in  ^  are  neglected.  This  model  is  based  on  the  notion  that 
turbulence  relaxes  to  Gaussian  statistics  for  the  large  scales  if 
inhomogeneities  are  removed.  It  is  therefore  better  suited  for 


conditioned  than  for 


Urt conditioned  moments.  The  pressure-rate  of 


strain  correlation  in  (13)  are  modelled  according  to  Rotta  [15]  ("return 
to  isotropy")  and  Hanjalic  and  Launder  [16] 

J  **  fZZ ) 


3 


with 


C,  -  /  S' 

cz  *04 


The  closure  of  the  interface  terms  S  is  based  on  the  following 

*7= 

considerations.  The  turbulent  zone  propagates  into  the  nonturbulent  zone 
by  viscous  transport  of  vorticity  into  irrotational  parcels  of  fluid. 
This  propagation  is  only  possible  if  the  net  effect  on  the  turbulent  zone 
stress  is  gain  at  the  expense  of  nonturbulent  zone  fluctuations  and  if  the 
fluctuations  in  the  turbulent  zone  are  weaker  than  in  the  nonturbulent 
zone.  Thus  the  model 


^  *  csr('-r>(w 


-  v*ii  ) 

*  ^  y  Ar 


SZS) 


emerges  with  -7  •  Thus  the  closure  of  the  turbulent  zone  stress 

equations  is  concluded.  The  constants  for  the  classical  terms  ^SC)C/tCz) 

are  taken  from  the  respective  references  and  the  constant  for  the 

interface  group  S*  was  established  by  computer  optimization. 

7® 


4.2.  Closure  of  the  nonturbulent-zone-stress  equations 

The  closure  for  diffusive  flux  and  pressure  correlations  for  the 
nonturbulent  zone  cannot  be  simply  carried  over  from  the  unconditional 
case,  because  the  nonturbulent  zone  fluctuations  behave  differently  (they 
are  irrotational)  in  the  limit  ^  >  C*  from  their  turbulent  zone 
counterparts.  Since  no  homogeneous  distribution  exists  in  the 
nonturbulent  zone,  the  decay  properties  of  correlations  in  this  zone  with 
distance  from  the  turbulent  zone  as  presented  in  (9)  are  used  to  estimate 
the  relative  order  of  magnitude  of  the  terms  in  (14).  Analysis  of  the 
stress  equations  shows  that  for  ^  O  diffusion  and  pressure  strain 


rate  correlation  become  the  leading  terms.  Furthermore  is  for  this  limit 


the  normal  stress  component  ZJ2X  dominant.  In  order  to  take  the  decay 
of  the  fluctuations  with  distance  from  the  turbulent  zone  into  account  the 
following  closure  expressions  is  used  with  a  time  scale  ^>7-  for  the 
non turbulent  zone 


) 

&*2  J 


and  with  the  decay  law  [17], [9] 


y 


-4 


where  y  denotes  the  normal  distance  from  the  center  of  the  turbulent 
zone,  we  obtain  for 


V 


where  denotes  the  turbulent  zone  scale,  the  decay  law 


/ 


3 


Let  furthermore 


be  approximated  by 


for  large  y 


,  then  follows 


for  rt  =  DC  *»3.  This  consideration  indicates  the  method  of  closure  for  the 
nonturbulent  zone.  The  classical  terms  will  expressed  in  terms  of  the 
unconditional  closure  with  the  modified  time  scale  2 «r  given  by  ( 26 ) . 


If 


The  dissipation  in  (14)  is  negligible,  because  viscosity  has  no  effect  on 
momentum  transport  in  irrotational  flow.  The  closure  for  the  diffusive 
flux  is  given  by  [12] 


9°  =  jf  7j[$  *  J  $ )]  "7> 

L  ?CL-*-S  *A  c  vi  /&  *  J 


Where 
A- 


4,„=  wsAqw  y  m) 


and 


A  -  s£  * 


(15) 


and  has  the  same  value  as  for  the  turbulent  zone.  The 


pressure-strain  rate  correlations  are  again  modelled  as  in  the 


unconditioned  case  with  the  modified  time  scale 


730) 


-Hr)  -*  +&)  *  4^  -H*)] 


The  definitions  of  7^  *2^  are  as  *n  (2^)  an<1  <24)  except 

nonturbulent  zone  quantities  and  averages  replace  the  turbulent  zone 
symbols.  The  second  part  in  (30)  representing  the  "rapid"  terms  is 
multiplied  with  the  intermittency  factor  to  modify  the  time  scale 
determined  by  the  strain  rate.  The  first  part  in  (30)  has  an  additional 
variation  with  in  terms  of  the  factor  77  ~f)  in  order  to 

compensate  for  the  decrease  with  of  the  second  ("rapid")  part.  Both 


r 
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constants  Ct  and  have  the  same  values  as  for  the  turbulent  zone 

counterparts  (22)-(24). 

_  0 

The  closure  of  the  interface  group  is  guided  by  the  properties 

of  the  solution  induced  by  the  second  part  of 

and 

Both  lead  to  wavy  solutions  and  thus  the  model 

sJ  *  &  +  y 5/ 


is  suggested.  It  is  a  preliminary  form,  which  gives  well-behaved 
solutions,  but  different  forms  using  the  consistency  condition  (17)  are 
currently  being  investigated.  Thus  the  closure  of  the  nonturbulent  zone 
stress  equations  is  concluded  without  new  constants. 


4.3.  Dissipation  rate 

The  dissipation  rate  £■  in  the  turbulent  zone  satisfies  a  complex 
equation  [16]  with  all  sink/source  terms  and  turbulent  diffusion  in 
non-closed  form.  The  closed  form  follows  Hanjalic  and  Launder  [16]  except 
for  the  turbulent  diffusion  Which  is  taken  from  Lumley  [12]  and  the 


interface  group 


with  C^m  0-7  determined  by  computer  optimization.  Thus  the  following 
equation  is  obtained: 


si 

9t 


j.  9 


fa  vje')  -  q, 


-is*g  *  c^/-r)f*^-Z)  nz) 

where  [ 12 ] 

~  ce  W 

The  constants  are  * /.  -2.-0  ,  •  <2 /S'. 


(5)  numerical  Solution 

The  system  of  partial  differential  equations  constituting  the 
conditional  closure  model  reduce  to  parabolic  form  for  boundary- layer- type 
flows.  Thus  they  can  be  solved  in  marching  type  integration.  The  present 
method  is  the  standard  finite-difference  procedure  developed  by  Patankar 
and  Spalding  [18].  The  new  aspect  of  the  solution  method  is  the 
introduction  of  a  block— solver  [19]  for  the  turbulent  and  nonturbulent 
zone  stress  tensors.  The  reason  for  this  is  the  strong  coupling  of  the 
stresses  as  a  consequence  of  the  diffusion  models  (19)  and  (27)  which 
makes  sequential  solution  likely  to  be  unstable.  Thus  are  the  stress 
equations  (13)  and  (14)  respectively  solved  simultaneously  forming  two 
block-tri-diagonal  systems  of  equations  with  blocksize  *  4  in  both 
cases.  All  other  equations  are  solved  sequentially  in  the  marching 


integration  step. 
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The  discretization  of  the  differential  equation  was  performed  using 
staggered  grids.  All  first  order  moments  ( intermit tency  factor  and  mean 
velocities)  were  defined  at  node  points.  Whereas  all  second  order  moments 
were  defined  at  the  midpoints.  This  leads  to  nearly  second  order  accurate 
representation  of  production  terms  in  the  stress  equations  and  the 
diffusive  term  in  the  mean  velocity  equations.  Furthermore  is  the 
stability  of  the  finite-difference  scheme  improved  [7]. 

The  initial  conditions  for  conditional  closure  models  can  be  set  in 
several  ways.  The  turbulent  zone  variables  and  the  intermittency  factor 
can  be  prescribed  corresponding  to  fully  turbulent  flow  or  the 
intermittency  factor  is  started  with  a  small  value  and  the  nonturbulent 
variables  with  laminar  profiles  representing  a  slightly  disturbed  laminar 
flow.  In  the  present  case  the  former  method  was  chosen. 

The  boundary  conditions  for  the  apparent  stress  tensor  in  the 
nonturbulent  zone  require  some  consideration.  Since  the  nonturbulent  zone 
stresses  cannot  approach  a  homogeneous  distribution  the  decay  laws 
discussed  in  ref.  [9]  are  used  to  establish  at  the  free 
boundary.  This  is  done  in  terms  of  a  gradient  condition 

37  f°'  y»yo 

Where  y  is  the  normal  distance  from  the  center  of  the  turbulent  region 
(symmetry  axis  on  location  of  maximal  shear  stress  in  turbulent  zone). 
The  location  of  the  origin  yo  for  y  is  not  known  exactly,  but  for  large 
values  of  y  this  relation  becomes  reasonably  accurate. 

All  calculations  were  performed  with  A/=  S~0  grid  points  over  the 
cross-section  and  the  number  of  steps  in  x-direction  ranged  from  800  to 
2000  depending  on  the  length  of  the  computational  domain. 


(6)  Applications 

The  conditional  second  order  closure  developed  above  was  applied  to 
the  calculation  of  several  plane  shear  layers.  All  calculations  were 
carried  out  with  the  same  set  of  constants  given  in  chapters  2.  to  4. 

6.1.  Plane  jet 

The  results  for  the  plane  jet  are  shorn  in  fig.  1.  to  fig.  10.  in  the 
nearly  self-similar  region.  The  experiments  of  Gutmark  and  Wygnanski  [20] 
and  Sunyach  [21]  are  used  for  comparison  with  the  calculations.  The 
intermittency  factor  profile  in  fig.  1  lies  between  the  two  sets  of 
measurements  (open  symbols:  [20],  full  symbols  [21]),  but  is  somewhat 
steeper  than  the  experimental  profile.  The  mean  velocity  for  the 
turbulent  zone,  the  nonturbulent  zone  and  unconditioned  are  compared  with 
the  experiment  [20]  in  fig.  2.  The  nonturbulent  zone  mean  velocity  o  in 
fig.  2  and  in  fig.  3  is  higher  than  the  turbulent  zone  mean  at  the  axis, 
because  the  nonturbulent  zone  shear  stress  (fig.  4)  is  always  less  than 
its  turbulent  zone  counterpart  thus  leading  to  slower  decay  of  U  The 
calculated  shear  stress  (conditioned  and  unconditioned)  in  fig.  4  is  lower 
than  the  experiment  [20]  in  the  outer  part  of  the  flow,  but  agrees 
reasonably  well  with  the  measurements  in  the  main  part  of  the  flow  field. 
Consequently  is  the  spreading  rate  dya f/ctjC  m  &  tOff  close  to  the 
experimental  value  €^£.5-  /a/x  >  O  //  .  The  normal  stress  components  are 
shown  in  fig.  5  to  fig.  10.  Their  relative  magnitude  can  be  evaluated 
from  fig.  S  for  the  unconditioned  case,  from  fig.  6  for  the  nonturbulent 


zone  and  from  fig.  7  for  the  turbulent  zone. 


It  is  clear 


I 

from  fig-  6  that  Vfz  becomes  dominant  as  y  increases,  which  is  in 

i  * 

I  accordance  with  the  relation  [9] 

i 

I  -  4^  +  &S1 

as  ^  — *  O  •  This  property  is  due  to  Lumley’s  diffusion  closure  (27)  and 
|  the  interface  terms  (31).  If  the  simpler  gradient-flux  model  of  Daly  and 

Harlow  [22]  is  applied  this  relation  cannot  be  satisfied  as  the  outer  edge 
is  approached.  The  comparison  of  the  normal  stress  profiles  with  the 
|  available  experimental  data  [20]  shows  reasonable  agreement  fig.  8  - 

fig.  10.  In  particular  are  the  shapes  of  and  well 

predicted. 

I 

6.2.  Plane  wake 

The  downstream  region  ®  200)  for  the  plane  waVe  of  a  cylinder  is 

compared  with  the  measurements  of  Fabris  [23]  and  Thomas  [26]  in  fig.  11  - 

1 

20.  The  profiles  for  the  intermittency  factor  in  fig.  11  is  about  ten 

percent  less  wide  than  in  the  experiments  (open  symbols:  [23],  full 

symbols  [24]),  but  the  slope  agrees  well  with  the  results  from  [23].  The 

mean  velocities  in  fig.  12  are  in  close  agreement  with  the  experiments 

Xj 

[23].  The  turbulent  zone  mean  velocity  U  shows  a  slight  bulge  near 
the  outer  edge.  Which  is  due  to  the  boundary  condition  set  to  be  equal  to 
the  free  stream.  The  nonturbulent  zone  mean  U  in  fig.  13  does  not 
allow  complete  comparison  because  only  limited  data  are  available.  The 
comparison  of  shear  stresses  fig.  14  and  normal  stresses  fig.  15  to 

I 

fiS  20  shows  much  the  same  properties  as  for  the  jet.  It  is  noteworthy  \ 

_  { 

that  tJy  and  1  in  fig.  19  are  close  to  the  measurements  as  for  ) 


the  jet. 
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6.3.  Plane  mixing  layer 

The  plane  mixing  layer  is  a  flow.  Which  is  rather  difficult  to 
predict,  because  of  its  sensitivity  towards  the  flow  conditions  and  the 
appearance  of  coherent  structures  [25],  which  is  reflected  in  a  wide  range 
of  observed  spreading  rates  and  significant  variation  in  the  experimental 
data  [26].  The  calculations  for  the  conditions  given  in  [27]  are  compared 
with  the  experiments  in  fig.  21  to  fig.  80.  The  intermittency  factor  in 
fig-  21  shows  a  broader  profile  than  the  experiments  of  Wygnanski  and 
Fiedler  [27]  but  agrees  on  the  high-speed  side  with  [28]  and  is  lower  than 
the  data  of  [29]  on  the  low-speed  side.  The  agreement  of  the  mean 
velocities  in  fig.  22  with  the  experiments  [27]  is  reasonable.  The 
turbulent  zone  profile  (full  circles)  from  [27]  is  not  approaching  the 
free  stream  on  the  low  speed  side.  It  is  not  clear  whether  this  is  due  to 
a  different  normalization  of  the  experimental  data  on  a  genuine  tendency. 
The  comparison  of  calculated  mean  velocity  for  the  nonturbulent  zone  with 
experiments  [27]  in  fig.  23  is  quite  good  and  extends  over  the  complete 
profile.  The  comparison  of  the  stress  components  in  fig.  24-30  with  the 
experiments  [27]  is  better  on  the  low-speed  side  than  on  the  high  speed 
side.  This  can  be  traced  back  to  the  dissipation  rate  equation  which  does 
not  produce  the  correct  length  scale  profile.  This  could  be  improved  by 
using  a  length  scale  which  is  constant  over  the  cross-section. 

(7)  Conclusions 

The  closure  of  the  first  and  second  order  moment  equations  for 
conditioned  variables  was  developed  and  the  resulting  system  of  model 
equations  was  applied  to  the  calculation  of  several  shear  flows.  The 
following  conclusions  can  be  drawn  from  the  properties  of  the  model  and 


the  comparison  with  available  experimental  data. 

(1)  The  intermittency  factor  plays  a  central  role  in 
conditional  closures,  because  all  unconditioned  moments  are  combinations 
of  conditional  moments  and  the  intermittency  factor.  The  source  term  of 
the  equation  for  the  intermittency  factor  requires  closure,  which  is 
constructed  as  difference  of  production  due  to  creation  of  apparent  stress 
in  the  turbulent  zone  and  the  inhomogeneity  of  the  g-  -distribution  and 
destruction  due  to  viscous  effects.  The  calculated  ^--profile  appears  in 
all  test  cases  in  reasonable  agreement  with  the  measurements. 

(2)  The  conditioned  moment  equations  contain  interface  terms,  that 
represent  the  transport  of  mass  and  momentum  due  to  the  movement  of  the 
interface.  If  it  is  assumed  that  the  conditioned  mean  pressure  in  both 
zones  is  equal,  then  follows  a  local  and  closed  consistency  condition 
eliminating  the  need  for  closure  of  the  interface  terms  in  one  of  the 
zones.  This  is  used  to  advantage  for  the  mean  velocity  equations,  but  for 
the  stress  equation  only  a  simplified  version  is  applied.  The  results 
show  that  the  nonturbulent  zone  mean  velocity  is  in  good  agreement  with 
the  limited  experimental  data  available.  The  comparison  of  the 
nonturbulent  zone  stresses  is  too  incomplete  to  draw  definite  conclusions 
however . 

(3)  The  notion  of  conditioning  allows  further  application  to  reacting 
flows  (in  particular  premixed  flames  [30])  and  generalization  to  nonlocal 
conditions  to  deal  with  structural  information  not  accessible  to 
single-point  theories. 
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Figure  captions 

Fig.  1  Intermittency  factor  for  the  plane  jet. 

<jcy/y0 y  Symbols:  0[20],  *[21]). 

Fig.  2  Mean  velocities  for  the  plane  jet. 

(Symbols:  [20]). 

Fig.  3  Mean  velocity  in  nonturbulent  zone  of  the  plane  jet. 

(Symbols:  [20]). 

Fig.  4  Shear  stresses  for  the  plane  jet.  (Symbols:  [20]). 

Fig.  5  Unconditioned  normal  stresses  for  the  plane  jet. 

Fig.  6  Nonturbulent  zone  normal  stresses  for  the  plane  jet. 

Fig.  7  Turbulent  zone  normal  stresses  for  the  plane  jet. 

Fig.  8  Normal  stresses  in  the  plane  jet  compared  with  experiment  [20]. 

Fig.  9  Normal  stresses  in  the  plane  jet  compared  with  experiment  [20]. 

Fig.  10  Normal  stresses  in  the  plane  jet  compared  with  experiment  [20]. 

Fig.  11  Intermittency  factor  for  the  plane  wake 

(  y^y/yos-  y  Symbols:  0[23],  »[24]). 

Fig.  12  Mean  velocities  for  the  plane  wake  (Symbols:  [23]). 

Fig.  13  Mean  velocity  in  nonturbulent  zone  of  the  plane  wake 
(Symbols:  [23]). 

Fig.  14  Shear  stresses  for  the  plane  wake  (Symbols:  [23]). 

Fig.  15  Unconditioned  normal  stresses  for  the  plane  wake. 

Fig.  16  Nonturbulent  zone  normal  stresses  for  the  plane  wake. 

Fig.  17  Turbulent  zone  normal  stresses  for  the  plane  wake. 

Fig.  18  Normal  stresses  in  the  plane  wake  compared  with  experiment  [23] 

Fig.  19  Normal  stresses  in  the  plane  wake  compared  with  experiment  [23] 

Fig.  20  Normal  stresses  in  the  plane  wake  compared  with  experiment  [23] 


Pig.  21  Intermittency  factor  for  the  plane  mixing  layer 

<  y*fy-yo*y/y„-y,.9) symbols:  0127],  •usj,  a  129] >. 

Fig.  22  Mean  velocities  for  the  plane  mixing  layer  (Symbols:  1 2 7  ]  > . 
Pig.  23  Mean  velocity  in  nonturbulent  zone  for  the  plane  mixing  layer 
(Symbols:  127]). 

Fig.  24  Shear  stresses  for  the  plane  mixing  layer  (Symbols:  [27]). 
Fig.  25  Unconditioned  normal  stresses  for  the  plane  mixing  layer. 

Fig.  26  Nonturbulent  zone  normal  stresses  for  the  plane  mixing  layer. 
Fig.  27  Turbulent  zone  normal  stresses  for  the  plane  mixing  layer. 
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Multi-scale  closure  tor  turbutent  shear  flows 
S.  Byggstoyl  and  w.  Kollmann 

1.  I ntroduction 

Turbulence  models  based  on  a  single  time  and  length  scale  have 
successfully  been  used  to  calculate  a  number  of  equilibrium  turbulent 
shear  flows  [1].  However,  when  the  turbulence  is  not  in  equilibrium  it  is 
generally  believed  that  it  is  necessary  to  include  more  than  one  length 
and  time  scale  in  the  modelling.  A  multi-scale  model  was  first  developed 
by  Hanjalic  and  Launder  [2].  Based  on  the  spectral  energy  transfer  they 
intuitively  constructed  transport  equations  for  different  regions  of  the 
energy  spectrum. 

In  this  paper  a  new  multiscale  concept  for  modelling  of  turbulent 
flows  is  developed.  The  concept  is  based  on  statistics  conditioned  upon  a 
positive  scalar  being  inside  some  given  intervals.  Each  interval  (zone) 
has  its  own  time  and  length  scale  and  different  statistics  can  be  used  for 
modelling  the  exact  equations  inside  each  zone.  The  dissipation  of  kinetic 
energy  is  taken  as  the  scalar  and  then  the  dissipation  term  in  the  equation 
for  turbulent  kinetic  energy  in  the  zones  is  viewed  as  an  independent 
scalar  variable.  The  modelled  transport  equations  for  each  zone  have  to 
be  complemented  with  either  the  equation  for  the  probability  density 
function  tor  the  dissipation  or  the  form  of  the  Pdf  has  to  be  given 
explicitly.  The  paper  is  organized  as  follows:  first  the  exact  transport 
equations  are  derived  and  discussed,  then  the  exact  equations  are  modelled 
and  solved  for  the  case  of  isotropic  turbulence  and  for  a  plane  jet  flow. 
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2.  DERIVATION  OF  EXACT  EQUATIONS 

First  the  exact  equations  serving  as  basis  for  the  Multi-scale  node! 
are  established.  These  Include  the  transport  equations  for  Interval ttency 
factor  and  conditioned  first  and  second  order  moments.  Some  of  the 
equations  can  be  found  In  Ref.  [3-6]. 

2.1  Intermlttencv  factor 

let  fry, l)  be  a  fluctuating  non-negative  scalar  satisfying  the  equation 

•il/'ri#)  /  Sf  <2-0 

where  S ^  expresses  the  production/destruction  of  the  scalar  (f>  (  5m  Is 
assumed  to  be  local  In  space  and  time  but  may  be  non-linear). 

Consider  a  discretization  of  the  range 

°  <  r,  <  <fz  ■  ■  ■  <</>»<«> 

and  let  the  Instantaneous  flow  field  <0  be  divided  Into  overlapping  zones 
)  ,M  Q  />  defined  as 


-9,  S  f£  >  <p/J 
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It  is  then  clear  that  Q-  ^  8  for  all  /  .  In  order  to  derive  exact 

t'ri  / 

transport  equations  for  the  flow  variables  Inside  each  zone,  a  sequence  of 
Indicator  functions  corresponding  to  the  sequence  Z^mJ  Is  defined 


as 


/: 


&(*>*)  >(f>* 

e>M  er+j/se. 


ZZ  2) 


The  surface  corresponding  to  each  in  the  range  &<f>)  is 

given  by  the  equation 


O 


Along  this  surface  dS^O  and  the  dynamics  of  the  surface  Is  described 
by  the  following  equation 

Q6Sa  +  -O  VS) 

where  ££  *  Is  the  instantaneous  velocity  of  the  surface  given  by 
Uc S‘ *  *  ^ 

where  A^*is  the  position  vector  of  the  surface. 

The  velocity  of  the  surface  is  further  split  into  two  parts 

Zl  ■=  +  fifr*  0i) 

where  Is  the  fluid  velocity  at  the  surface  and  is*  is  the 

propagation  velocity  of  the  surface  relative  to  the  fluid.  *1*  is  the  unit 
normal  vector  on  S*  defined  as  positive  pointing  into  the  domain  2>K  . 

The  indicator  function  is  obviously  constant  along 

>^1  Z )  and  satisfies  the  equation 
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l'.v- 

EvTj ■ 


3tIA  /  -O 

Using  (2.4)  this  tranforms  Into 


ft.  5) 


ZT<  +  V*  Z1* 


m.  Je 


f2  t>) 


The  derivatives  of  the  indicator  function  can  be  expressed  as 


ft  7) 


flj;  -  is  S& 

*f  *  *  K  OO 


■>P2<  Zs^rs*) 

where  <T<  )  is  the  Dirac  pseudo-function.  Introducing  (2.7)  in  (2.6)  gives 


By  taking  the  mean  value  of  this  equation  results  in  the  equation  for  the 
intermittency  factor  3  ( } 

A 


A  * 

=  <'y*/P5t/Jrs/cj> 

ft.  a) 

s  point  conditional  mean 

and  fluctuating  vel 

ocities  are 

introduced  as 

^  u 

Zj  * 

*  * 

~JL, 

-  ZJ* 

$C 

A 

)  * 

0C 

fZ  3) 

<0-1.)!+) 

0  > 

v  *  V 

)  K 

SHjl. 

where  D*  is  the  mean  velocity,  at  )  of  fluid  belonging  to  SL  and 


ZjJ1  is  the  mean  velocity  of  the  fluid  which  does  not  belong  to 


Am  l  - - >  J 
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.  . . .  I  |  |  |  nil  I  I  I  mi . ,  ■'■■■!  I.  _ 

By  Introducing  this  Into  (2.8)  gives  the  following  equation  for  the 
intermlttency  factor 

■+  %<'%*#■.)  ‘  <*■»> 

The  relative  propagation  velocity  V  can  be  expressed  In  terms  of  the 
discriminating  scalar^  by  using  (2.3)  and  (2.4) 

4 9s  +  ZfVc?)  -  f'Vr#/ 

which  Is  valid  at  the  surface 

By  Introducing  (2.1)  )/*  can  be  expressed  as 

fuki  ^  )J^o  "■  ^ 

which  shows  that  i/K  is  not  bounded  for  /Pr0/£^  — ?  O  ■ 

By  introducing  this  into  (2.10)  the  intermlttency  equation  takes  the 
following  form 

~  <'f*LfrQc4>)+  (2.V) 

The  right  hand  side  of  this  equation  will  be  denoted  by  and 
expresses  the  mass  entrainment  per  unit  mass  of  fluid  from  zone  k~!  into 
zone  k  .  It  can  be  shown  that  this  term  is  positive  or  negative  depending 
on  the  nature  of  the  scalar  and  the  threshold  value  (f>K  . 

In  order  to  develop  equations  for  the  flow  variables  inside  zones 
where  the  value  of  $>  is  inside  some  specified  interval  a  non-overlapping 
decomposition  must  be  used.  This  Is  done  by  defining  2/it)  as 

~  Q  ~  2^ 


which  satisfy  the  requirement 


%{*)  =  $t>  for  k  4  L 


The  appropriate  Indicator  function  and  Intermlttency  factor  for  this 


decomposition  re,  respectively 


“  j\c  ft*/ 

The  conditional  mean  and  fluctuating  velocity  for  fluid  belonging  to  bfte) 


/t'3) 


Is  then 


/*)  j 


>  %*•  K  -  v™ 


By  using  (2.14)  and  (2.12)  gives  the  following  equation  for 


5  '*)  St'S) 

Note  that  the  divergence  of  £j^1s  non  zero  and  the  relations  f efts)  * /FS/J(S)  ) 


By  multiplying  the  Instantaneous  momentum  equation 

with  the  Indicator  function  X(cC*'}4)  and  averaging  gives,  after  some 
operations  with  the  Dirac-pseudofunction  ,  an  equation  for  the  zone 

z’jf) 

conditioned  mean  velocity  ZJL  . 


-  v*  -fi„><$r4r>J  *  *€-  *C' 


The  momentum  fluxes  through  the  iso-scalar  surfaces  and  SM4t  are 
given  by 

>£&,/»  *  i/*~  v&  H  '£+  §n^)^sK)>  -  g  wn 

f2  *{,) 

2,  **• 

and  follows  by  increasing  the  index  k.  The  intermittency  equation 

(2.  IS)  allows  rearrangement 


f* 


a.  17) 


Comparison  of  (2.17)  with  the  equation  for  the  unconditioned  mean  velocity 

reveals  the  appearance  of  two  new  term  groups.  The  first  group 
proportional  to  the  zone-conditioned  Reynolds-stress  arises  in  the  elimination 
of  the  intermittency  factor  on  the  left  hand  side.  Its  role  can  be  elucidated 
for  thin  shear  layers  like  jets,  mixing  layers  and  boundary  layers  ano  low 
values  of  k  (say  one)  and  O  .  Then  can  be  seen  that  this  source 

pushes  the  conditioned  velocity  profile  further  out  than  the  unconditioned, 
hence  increases  the  spreading  rate.  The  second  group  is  a  collection  of  point 

statistical  moments  representing  momentum  transfer  through  the  interface. 
Chapter  three  is  devoted  to  the  properties  of  this  group. 
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2.3  Zone-averaged  Reynolds- stress  tensor 


The  transport  equations  for  higher  moments  of  zone-conditioned 


quantities  can  be  obtained  without  difficulty  using  the  properties  of  the 


indicator  functions.  The  Reynolds-stress  tensor  satisfies  for 


zone  (ktk+!)  the  following  equation 


*  V* 


f4c)  -  0 


+f  44~j)  -  Sv  sUa,HfA,> 


/2J8) 


The  Reynolds-stress  is  defined  by 


and  the  turbulent  flux  of  by 

7® 


Conditional  statistics  introduce  two  new  term  groups  as  in  the  mean 


velocity  equation  (2.17).  The  first  group,  arising  from  the  elimination 


of  the  intermittency  factors  on  the  left  hand  side,  can  be  recast  as 


+  fr-XL.  >>J  -  fa  -fi*)-  Kf^fa  -A- ) 


+  viscous  +  Divergence  terms. 


J>'*>  A  ,*,') 

where  a  %  +  4-  ‘ 

This  term  group  can  be  expected  to  provide  additional  transport  of 
thus  changing  spreading  rates  of  -profiles  in  thin  shear 

layers.  A  more  detailed  estimation  of  its  effect  is  however  difficult  due 
to  the  lack  of  experimental  information  on  the  i ntermittency  sources.  The 
second  group  of  point-statistical  correlations  defined  by 

s  +J-  <p fa )> 


will  be  discussed  in  chapter  three. 

fM) 

The  equation  for  the  kinetic  energy  of  turbulence  k  for  the 
-zone  deserves  special  attention  for  £  £  .  It  follows  from 
(2.18)  as 

-  e  *> 


+  +  v^k't:)  ^A^fa  "A"  ) 


fi-  A 


[Z  -f*„  )  -  A'Vfe* »  r  (^u,rfsw  )>J 


-*T‘ 


The  dissipation  rate  conditioned  for  zone  k  is  obviously  bounded  by  Q 
and  (pkfl  ,  hence 


10 


£'*’  -  jf9<  +  )  +  Ofa  -  <?„,) 

If  the  interval  length  —  5^,  /  is  small  enough,  this  does  not 
constitute  an  unknown  correlation  (for  which  an  equation  could  be  derived) 
but  an  additional  discrete  independent  variable.  The  zone-conditioned 
moments  can  be  considered  therefore  functions  of  and  (^>  .  If  we 

assume  for  the  moment  that  r~.u  >0  and  r~  >  O  we  see  that 

«*  toe  ^ 

L 

several  inverse  time  scales  can  be  formed  for  the  /r^zone: 


Energy  input  from  zone  k-/ 


k'M)  +k'*") 


F 

<x*c 


"  }* 


HH 


Energy  output  to  zone  k+) 


£ 


-/ 


F 


*+/ 


k'k)+  k 


fk+n 


<}*- 


Energy  transfer  through  zone  k 


/  /F~  F"/ 

fa- 


Average  scale  for  zone  k 


v-/ 

ZkfM) 


The  input  and  output  scales  make  sense  however  only  for  finite  intervals 
Afa.  Hence  only  the  transfer  scale  and  the  average  scale  will  be  used 
in  the  following  together  with  the  corresponding  length  scales  formed  with 
appropriate  powers  of  k  . 
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2.4  Relation  to  unconditional  moments 

Unconditional  moments  can  be  recovered  from  the  sequence  of 
conditional  moments  via  local  relations  following  from  the  definitions 
given  at  the  beginning  of  this  chapter.  Noting  that  m  1  and  u 

we  obtain  for  the  unconditional  mean  value  of  a  fluctuating  quantity 

AJ 

ip  -  Z  -  fa, )  t:  w) 

k»o 

for  the  special  case  A/*/  we  note  that  /V  and  and  then 

reduces  (2.2U)  to  the  well-known  relation 

Z  -  +  f'-r*? 

Furthermore  is 

Z  (fa  -fa*,)  -  1 

Jc~o 

and  the  terms  in  the  sum  are  nonnegative  due  to  monotonicity  of  the 
intermittency  factors.  For  second  order  moments  follows  similarly 

wz  =  z  <0n’Z">z*  -  >  +z  zyz?*  -/**>) 

&=o  i<x° 

-zz  zyzp*  y  -fa,  >  f2 2,) 

k*-0  L-O  v 

For  the  special  case  N  =  1  it  reduces  to  the  relation 

ZZ’  '  fa  +  pif/-pW'>-,PmPy’",-y>'a) 

The  formulae  for  higher  order  correlations  can  be  derived  without 
difficulty. 


Zone-cond i t i oned  moment  equations  contai''  a  tenr  group  amount 


the  transport  effects  through  the  iso-  sea  ar  surface  represented  tv  tnu 
impl icit  equation 


-  fit*,**  -  </>«  -o 


This  tenn  grojp  appears  in  the  intermittency  tdutor  equations  as 

IM°  =■  ))  rs-D 


and  in  the  equation  for  moments  of  order  one  and  higher  as 


jj  *  _  U**1 
X*  3  — 

*  /* 


n  2) 


provided  the  demon inator  is  nonzero.  denotes  a  collection  of 

point-statistical  correlations  of  the  arguments  velocity  and  pressure  with 
the  relative  progressive  velocity  V  and  the  normal  vector  of  tne 

iso-scalar  surface  S*0. 


3.1  Source  terms  of  intermittency 

The  term  group  (3.1)  represents  the  transport  of  the  iso-scalar 
surface  itself  and  contains  productive  and  destructive  effects.  This  can 
be  seen  by  considering  the  limit  First  we  notice  that 

tO  60 

/*-/«/  *  ))  “f&y  ~  [dy  Pf?) 

% 
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where  Pf(f,rt4)  is  the  one-point  pdf  of  the  discriminating  scalar  0  . 
For4y2r>0  but  sufficiently  small  and  P(ip)  sufficiently  smooth  follows 

fa  -  fat  -  &</>*  Pty*)  +  hot*  {3.3) 

Furtnermore  follows  from  the  definition  [2/b) 

vft)  =  Su/<t  =  <A>  #.*; 


which  is  the  mean  of  V  subject  to  the  condition  that  the  scalar  ^assumes 
at  the  same  (fixed)  location  the  value  .  The  relations  (3.3)  and  (3.4) 
are  the  primary  tools  in  analyzing  the  sources  (3.1).  Making  the 
^-interval  sufficiently  small  and  applying  (3.3)  and  (3.4)  to  (2.14)  we 
obtain 

+  o 

Dividing  the  and  letting  *%  -*0  leads  to 


* 


(3  S') 


The  limit  on  the  right  hand  side  can  be  evaluated  by  direct  derivation 

A 

the  pdf-equation.  Defining  the  fine-grained  pdf  P  by  [7] 

p  =  ,  p  -  <p> 

the  pdf  equation  follows  [8]  as 

DtP  *  2 ft U/0-p>P)  •~£<'eaP>  -  v  r 

* 

Equating  the  right  hand  sides  of  (3.5)  ana  (3.6'  \ 

result 

7 


A***  A0  ]'b 
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Note  that  £.s  Is  defined  as 

es  m  r  v0.v#> 

From  (3.7)  follows  that  the  limit  on  the  left  hand  side  exists  and 


therefore 


2.1'yrtsj)  -  XCSrpP)  +  ZLse'J*} 


and  finally 

)  -  Sry)  Pf</)  +  (e  *P>  +  tons/  S3  ff) 

Equation  (3.8)  allows  a  detailed  discussion  of  the  interim ttency  source/ 
sinx  term  for  the  case  Pfif)  close  to  1 09- normal.  First  we  note  that 

sr'p  >  >  o 

and  therefore  will  change  sign  at  least  once  because 

~*C  as  ana  •  0  .  T he  second  derivative  in  (3.7) 

■  1  '•  be  negative  around  the  mean  {#}  of  ?  ,  Sj,  hence  is  isett  ing  Sm  O 
for  t  ne  moment  •  y  '>>0.  For  va’ues  tar  away  from  the  second 

Je r  1 » a  t  ■  ve  1  n  w  1  >  1  oe  pos  1 1  ’  ve  maa  1  ^ C  1  f  P< ''(f )  1  s 

smooth  •  •  t  n  a  v  •  nq  *  e  pea*  n»>r  ,  ;t  •  ?  specialise  the  d  1  sc  r  lmi  nat  1  nq 

a  4’  t.  (pm  £  ,  the  Jss'pation  'ate  of  t  .r;.u  ent  » '  net  1 .  energy,  and 


o'i  v  >e  ’  a  •  eg  ■  nr 


the  turbo  ent  f  om  f  t  id  wnere  ,  the 


t  t  ne  '  ■  cs  ipat  ■  or.  of  & 


a  v,/c 

\  *  e 


^  >  r 


V  •  y  if  '  iro 


,  a  ' 


c  e  Cf  ^  a 1  '  ’ T  a '  * 


*%*■*•;  \  v 


1  '  r  '  Vme, 


»  a  >  •  •  at  rme*  ?  -  «■  '  e  '  a,  » 


v  t  *  a  1 


’*  .r  \y hr  *r  ^ 
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for  (p*  and 


< v* >  >  o 


fory^yfS)'  The  scalar  dissipation  terms  appearing  in  (3.7)  and  (3.8) 
are  sketched  in  fig.  1  for  the  case  of  F*«p)  being  close  to  log-normal 
form. 

The  effect  of  source  Sf<f)  of  the  dissipation  rate  £,  on  the 
intermittency  factor  becomes  apparent  from  (3.8).  If  Sff&J  > O  then 
Y*  is  produced,  but  SapJ  JO  implies  at  the  same  time  a  shift  of 
towards  higher  values  of  (fi  and  the  growth  of  X*  is  changed  according  to 
the  shape  of  F*  . 

The  shape  of  the  ^ -prof  i  les  can  be  inferred  from  (3.5),  (3.8)  and 
the  boundary  conditions  with  respect  to  the  ^>-axis. 

J-C  -  >  for  ,  /*„  -  O  for  f  ■ 

Equation  (3.8)  shows  furthermore  that  (provided  the  constant  in  (3.8)  is 


zero)  the  real lzabi 1 ity  condition 


O  £  yk  <  /  for  , 


is  satisfied,  because  the  right  hand  side  is  the  source  term  of  the 
equation  for  (integrate  (3.6)) 


If)'!/) 


T 


which  satisfies  O  £  tj>  £  /  . 

For  small  <J>M  (close  to  the  boundary  <f0m  O  where  *  /  )  the 

profile  of  ^  will  be  close  to  ft0~  I  if  the  turbulent  Reynold s-number 
is  high,  because  the  pdf  of  £.  will  be  small  for  small  values  of  and 
be  concentrated  at  high  values  of  <j)  and  therefore  Sf<p)  and  will 

be  small  too.  For  high  values  (/>M  <  the  profile  of  ^  will  peak 

"ear  where  the  source  of  may  still  be  positive.  Further  away 

from  the  location  of  )  where  (£)  the  source  of  t/  will  be 

away  M. 

negative  (see  fig.  1).  The  sequence  of  intermittency  factors  is  however 

monotonical ly  decreasing  everywhere  in  X  -space. 


Figure  1  -  Distribution  function  and  scalar  dissipation  function 
for  two  different  spatial  locations  with 
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These  properties  of  are  summarized  In  fig.  2  for  the  example  of  a 
thin  shear  layer 


Figure  2  -  Intermittency  factors  for  plane  jet  and  N-4, 


It  should  be  noted  that  for  sufficiently  small  the  profile  of 

corresponds  to  the  profile  of  the  classical  intermittency  factor  which 


distinguishes  between  turbulent  and  non-turbulent  states  of  the  flow.  A 


large  amount  of  experimental  data  is  available  for  this  case  (references 


given  in  [6]) 


The  relative  progresson  velocity  r  of  the  iso-scalar  surface  S  m O 


was  derived  in  sec.  2.1  as 


V«f>,  -  +  Sr<p) ] 


fZ  3) 


<p.y 
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This  can  be  applied  to  (3.3)  and  leads  to  a  new  expression  of  the  scalar 
dissipation  term  of  the  pdf-equation  for  a  fluctuating  scalar,  we  find 

V 

<fe%p  >  *  -</*>> 

-•0 

Hence  Is  the  scalar  dissipation  of  P  at  the  integral  of  the  relative 
progression  speed  of  all  iso-scalar  surfaces  with^V^J.  Equation  0.9) 
shows  now  that  the  Instantaneous  values  of  i/  are  not  bounded  and  the 
amount  of  scalar  dissipation  of  P  at  depends  on  the  Pdf  of  (p  . 

<Cesfi )  -  Jo/(/>)  f  -^>>)}  (3*0) 


-•O 


3.2  Source  terms  for  zone-conditioned  mean  velocity 


The  term  group  (2.16)  represents  the  transfer  of  momentum  through  toe 
iso-scalar  surface  SM*  O  and  contains  correlations  of  zone-conditional 
velocity  with  the  progression  speed  of  the  surface,  the  zone-conditional 
pressure  with  the  fluctuating  normal  vector  of  the  interface  and  viscous 
interactions  of  velocity  and  normal  vector.  The  limit  -*0  revea  1  s 

again  the  effect  of  this  group  on  the  zone-conditional  moment  concerned. 
With  (3.3)  we  can  write 


<  « 


fi  o  ( 


The  terms  constituting  H  contain  a  viscous  contribution  of  the  form 


<3,  A'  <VK  n*  ) > ) 

which  is  negligible  compared  to  the  other  terms  for  high  Re^  -numbers 
because  it  scales  with  the  mean  field.  Hence  it  will  be  neglected  and 
appears  to  be 


\ 


*  Pf<p« ] <1 V//0 m<pH > 

T"  ^ 

Then  follows  for  JL K 

_*  P/</>« )<[  Jjy *</>«)  —  Ps<p+  * &<f* ) <f  JA 0 =9>«+W> 


«  Ay*  Pf(/>A  ) 

The  fluctuating  components  are  in  the  limit  &%-*o 

tj'k'  -  ZL-  <'VK/$ 


T*  * 

and  J 1  turns  out  as 


•C-  -Tsr, 

rs.  if) 

provided  Pf(f0t)>D.  The  terms  do  not  contribute  to  the  overall 

(unconditional)  momentum  balance.  This  can  be  seen  from  the  expectation 
Of  Xx 

•o 

<T*>  ■=  J 1?  Pr<f>  V-1 


With  (2.32)  we  obtain 


/"-A? 

=  /^y>;  //^  m(/y /  ~  & 

</>=o 

Hence  transfer  the  point-statistical  terms  -C*  momentum  along  the  y?-axis 
without  affecting  the  unconditional  momentum  balance.  Since  all 
zone-conditioned  mean  velocities  must  satisfy  the  same  boundary 
conditions,  momentum  is  transferred  from  the  boundaries  to  all  zones.  For 
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small  values  of  <p  the  zone-conditioned  shear  stress  can  be  expected  to  be 

small  and  will  therefore  be  insufficient  for  momentum  transport  to  keep 

the  corresponding  zone-mean  velocity  field  stable.  Hence  must  the  point- 

statistical  term  group  transfer  enough  momentum  out  of  these  zones  to 

ensure  stability.  As  increases  the  Reynolds-stress  in  these  zones  will 

'  —.K 

become  more  important  for  momentum  transfer  and  ^  will  decrease  and 

—  K 

change  sign.  For  high  values  p  the  role  of  JL ^  will  be  opposite  to  its 
effect  for  low  CP  . 


3.3  Source  terms  for  zone-conditioned  Reynolds-stresses 


The  term  group  (2.19)  represents  the  transfer  of  stress  across  the 
iso-  scalar  surface  O  .  Its  properties  are  similar  to  the 

corresponding  group  in  the  momentum  equations.  In  the  limit  O  we 

f  ind 


if  -  i- 

+  v“'«* )  -*>%  fv^upn^]/#  -  >/  f3- 


where  the  last  viscous  term  in  (12)  is  again  negligible.  The  unconditional 


iro 

t c 


expectation  of  V  is  zei 

/ )  =-  Jdy>  Pry)  ~  & 

o 

indicating  that  the  point-statistical  group  (2.19)  transfers  stress  and  in 
particular  kinetic  energy  along  the  <^-axis  without  affecting  the 
unconditional  balance.  Furthermore  must  change  sign  at  least  once 

with  respect  to  <f).  If  the  choice  <p  s  t  is  considered,  the  high 
values  of  £  correspond  to  the  small  scales  of  the  turbulent  motion. 


For  high  turbulent  Re-numbers  local  isotropy  can  be  assumed  for  this  range 
and  small  normal  stress  levels.  Then  follows  that  there  is  little 
production  of  zone-Reynolds-stress  due  to  mean  strain  rate  and  most 
production  will  be  supplied  by  the  transfer  term  (3.12),  which  is 
therefore  positive  for  high  (^-values  and  and  negative  for  small 

^-values. 

3.4  Multiscale  closure  models 

The  first  step  in  the  development  of  mutli-scale  models  based  on 
conditional  statistics  is  the  choice  of  the  discriminating  scalar  . 
The  most  natural  choice  is  the  dissipation  of  turbulent  kinetic  energy 
because  this  quantity  appears  in  the  equation  for  turbulent  kinetic  energy 
(and  in  the  equation  for  the  individual  normal  stresses).  If  the 
turbulence  field  is  close  to  equlibrium  it  is  probably  not  necessary  to 
solve  an  equation  for  the  PdF,  but  instead  use  a  prescribed  Pdf  expressed 
in  terms  of  some  of  its  moments.  If  £  is  taken  as  the  discriminating 
scalar  then  the  log-normal  distribution  may  be  a  good  choice  [9],  In  a 
non-equilibrium  situations  it  is  necessary  to  solve  a  modelled  equation 
for  the  intermittency  factors  fa  or  to  solve  the  Pdf  equation. 

4.  Closure  of  the  Multiscale  Equations 

In  this  chapter  closure  of  the  multiscale  equations  is  discussed  and  a 


simple  closure  model  is  suggested.  The  model  is  applied  to  decaying 
Isotropic  turbulence  and  to  a  plane  jet. 


4.1  Multi  scale  Model  for  Decaying  Isotropic  Turbulence 


In  the  case  of  isotropic  turbulence  the  equation  for  zone  conditioned 
turbulent  kinetic  energy  simplifies  to 


<2  k 
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8y  using  (3.8)  and  (3.12)  this  equation  transforms  into  (by  letting 

&(p*A €->Q. 

4 f Pee j)  .  -Qe?=  -  ar*'pfe)  r*-2) 


where  F  s  —  (See  (3.12)) 

Equation  (4.2)  has  the  same  form  as  the  equation  for  the  energy 
spectrum  for  isotropic  turbulence.  This  is  seen  by  first  integrating  over 
all  values  of  E  ,  then  the  following  equation  is  obtained. 


d(  <*>  -  -  <e> 


which  is  the  equation  for  kinetic  energy  for  decaying  isotropic  turbulence, 
tty  integrating  (4.2)  from  0  to  some  yields 

%k,  -  "  FfS, )  -  Fy  f«.3) 

where  k,  — jds  krB)  >  E,  *  E  P(E) 

o  c 

the  term  -  expresses  the  net  energy  transport  from  zones 

with  E<  Et  to  zones  with  E  >  E,  f  this  is  illustrated  in  Fig.  3. 
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Figure  3.  -  Schematic  illustration  of  the  energy  transport  (eq. (4.3) ) 

As  seen  from  (4.2)  the  Pdf  of£,  ^fag^is  needed.  This  can  be  obtained 

by  a  direct  solution  of  a  modelled  transport  equation  for  the  pdf  or  by 

prescribing  a  pdf  which  depends  on  some  of  its  moments.  Both  of  these 

methods  will  be  discussed  next. 

4.1.1  The  Transport  Equation  for 

The  exact  transport  equation  for  PfE.)  can  easily  be  derived  from 
the  Navier  Stokes  equation  and  its  form  is  given  by  (3.6).  The  equation 
for  is  derived  in  Appendix  1  and  for  the  case  of  isotropic 

turbulence  i±  is  given  by 

-  -  £ <”?£■?££>  -  te 

+  *  '*■*> 

where  £  is  the  probability  space  variable  corresponding  to  the  scalar 

£  *  *>  2*1 
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-*  i—  -  - 


to, 


The  scalar  dissipation  term  is  modelled  according  to  [8] 
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where  the  transition  probability  T~  is  given  as 
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and  the  timescale  2~  as 
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It  is  easily  shown  that  the  expression  (4.5)  decreases  the  variance  of 
the  Pdf  without  affecting  the  mean,  therefore  the  last  term  on  the  right 
hand  side  of  (4.4)  must  be  responsible  for  the  decay  offe).  This  term  is 
modelled  as 

-*v  %  fa  *  **£*&£%  = 

=  -p  OgfeP)  hs) 

which  represent  a  convection  of  the  Pdf  towards  B=0  ( for  Cr>0) 
with  convection  velocity  £^.5  ^  .  Because  this  velocity  is  not 

uniform  but  depends  on  £.  the  higher  order  moments  will  also  be  influenced 
by  this  term. 


The  modelled  transport  equation  for  P(£)  is  then 


<?  O 


By  multiplying  (4.9)  with  £  and  EZ  respectively  and  integrating 
over  £.  transport  equations  for  and  are  found.  They  are 

c/i<e>  * 

+2ez)<'£,l>  (*'») 

Equation  (4.10)  is  the  same  as  the  )  equation  in  the 

model  of  turbulence  (in  the  case  of  decaying  isotropic  turbulence)  and  the 
constant  Cz  is  therefore  put  equal  to  Czm1.92.  It  is  known  from 
experiments  that  the  ratio  <€*>/&>  increases  with  increasing 

Reynolds-number  of  turbulence,  then  by  solving  the  Sk}  ~  model 
together  with  (4.11)  it  can  be  shown  (See  Appendix  2)  that  this  constant  A 
can  be  given  when  comparing  (4.11)  with  Kolmogorov's  modified  equilibrium 
theory  [10]  to  be  discussed  next. 

4.1.2  Prescribed  form  of  the  Pdf 

In  order  to  take  into  account  the  intermittent  behaviour  of  the 
dissipation  Kolmogorov  [10]  and  Obukov  [11]  refined  the  "universal 
similarity  hypothesis"  [12]  by  postulating  that  the  dissipation  of  kinetic 
energy  was  log-normally  distributed  with  a  variance  dependent  on  the 
Reynolds  number.  More  specifically  they  postulated  that  the  Pdf  of  the 
dissipation  was  given  by 

Pl'e>  =  7*7 A  e  ^ ^ ] 

where  m  =»  }  AZ  m  £++  e) 
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The  parameters  Att  and  can  be  related  to  the  mean  and  variance 

of  £.  by  the  following  expression.  [See  Appendix  3.] 

« ' u '*,s) 

Further  postulated  (Kolmogorov  [10]  and  Obukov  [11])  that  the  variance 
of  &  was  locally  related  to  the  macro-  and  dissipation  length  scales  of 


the  turbulence  by  the  formula 


i/or  / (A* —  A  + 


/s) 


where  A  is  a  "constant"  determined  by  the  macro  structure  of  the  fio* 
and  is  an  "universal"  constant.  l~0  is  the  macro-length  scale 
determined  from  the  integral  of  the  two  point  correlation  function  ana  *j 
is  the  dissipation  length  scale  defined  as 
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introducing  (4.13)  in  (4.15)  gives  the  following  expression  r^r  tne 


variance 
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By  asssuming  that  «>  *  *)  ana  def  I  n  i  ng  a  t  Vj 

Reynolds  number  as  V)  (4.16)  can  be  rewritten  as 

=  <fe>7 e\  -/  ] 

Masiello  [13]  detennined  experimentally  the  constants  a  ana  t 
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From  the  work  of  Driscoll  [14]  the  constant  OC  is  estimated  as 
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4..  Agp  utKT  of  mu  ti^ale  modelling  to  turbulent  shear  flows 

*ne  qeret a  equation  for  /one  conditioned  turbulent  kinetic  energy 
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«*nere  Fl.  .  is  the  ncn-viscous  part  of  W~  A  .  When  solving  an  equation  for 


*4“  the  only  new  terms  to  be  modelled  are  the  production  teri.i 
'and  the  turbulent  transport  term  .  However,  in  this 

mitial  staqe  of  model  development  the  zone  conditioned  mean  velocities 


1  oe  assjmed  equal  to  the  uncond 1 1 lona  1  mean  velocity. 
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T ne  zone  conditioned  shear  stress  is  assumed  to  be  proportional  to 

tne  urtcona i t loned  kinetic  energy  and  the  following  closure  model  is 
suggested . 

t'A ’  -  C  -%£-  *  *<%>)  -  i&*>  ft-Oi 


The  transport  term  [P  1  s  modelled  according  to  the  turbulent  viscosity 


concept 


(4,21) 


The  modelled  equation  for  zone  conditional  kinetic  energy  for  the  case  of 
a  two-dimensional  turbulent  shear  flow  then  becomes 

l/kMP)  +fu>0jJcmP)+  <V> 9y "P)  fa ^P)] 

-  e'uP  +  c,,  +  c„('etv-f£>)P>  f*.&) 

From  this  equation,  two  zones  are  constructed  by  Integrating  over  £  from 
®nd  from  €t-+»o .  The  kinetic  energy  in  the  two  zones  are  defined  as 
<£*/  *> 

•  Jo/e  k  ,  k£  * Jvte  k'Vej  Fte) 

O  €* 

V  ^ 
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In  cases  where  the  mean  dissipation  SS)  has  a  large  spatially  variation 
it  may  not  be  possible  to  find  a  representative  value  of  £,  for  the  whole 
flow  field.  Then  it  is  better  to  express  £t  as  a  fraction  of  <£)  and 
transform  tne  equations  (4.25)  and  (4.26)  accordingly.  This  transformation 
is  shown  in  Appendix  4. 


5.  Results 


5.1  Decaying  isotropic  turbulence 


In  order  to  solve  eq.  (4.24)  the  following  dimensionless  variables  are 
introduced 
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where  and  the  subscript  o  indicates  the  value  of  the  variable 

°  <€>o 

when  imO. 

Equation  (4.24)  transforms  then  into 
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The  initial  Pdf  is  taken  as  locj-normal  and  the  distribution  of  k  as  a 
function  of  f  is  assumed  proportional  to  r^f£)  . 
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which  satisfies  the  requirement 
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Equation  (5.1)  is  solved  by  using  both  a  prescribed  pdf  and  by  solving 
the  pdf-transport  equation.  For  the  results  shown  in  figs.  5-8  the  pdf  is 
prescribed  as  log-normal  and  Kolmogorov's  expression  for  the  variance  is 
used.  The  mean  dissipation,  which  is  a  parameter  of  the  pdf,  is 
calculated  from  the  standard  fit )  -(s)  model. 

Fig.  5  shows  the  distribution  of  as  a  function  of  £  for  three 

✓  ✓  w- 

different  times  £  -  0,  0.21,  0.42.  Fig.  6  shows  the  "spectrum"  , 

fig.  7  shows  the  interzonal  transport  of  energy  F  and  fig.  8  shows  the 
evolution  of  the  pdf.  The  results  are  in  qualitative  agreement  with  what 
one  would  expect,  but  due  to  the  complete  lack  of  experimental  data  it  is 
not  possible  to  draw  any  conclusion  about  the  quantitative  behavior. 

Figs.  9-12  show  the  corresponding  results  when  eq.  (5.1)  is  solved 

together  with  the  transport  equation  for  PfC )  .  The  constant  A  in  the 
pdf  equation  is  calculated  from  eq.  (8)  in  Appendix  2  and  Ct  is  put  equal 
to  Ct*  O.  /  .  The  results  from  this  calculation  are  nearly  identical  witn 
the  results  shown  in  figs.  5-8. 

5.2  Plane  jet 

Tne  results  from  the  calculation  of  a  plane  jet  are  shown  in  figs. 
13-16.  The  value  of  the  discriminating  scalar  C ,  that  defmesthe  two 
zones  was  set  equal  to  £  *  10*"^  m^/s^.  The  exit  velocity  and  tne 
nozzle  width  were  35”  m/  s  ana  J)  *  3  x  lO~m,  respectively.  Figures 

*Y/  ^ 

13a-c  show  a  typical  development  of  the  zonal  kinetic  energies  A,  and  A^  . 

Figure  13a  shows  that  for  x/D  »  10  the  main  part  of  the  kinetic  energy  is 

contained  in  zone  2  (  £>E,  ).  Further  downstream  the  dissipation  rate 

decreases  and  more  and  more  of  the  energy  is  contained  in  zone  1  (£<■£,). 

For  X/D  *  60  fig.  13c  zone  2  is  drained  out  and  all  the  energy  is 
contained  in  zone  1.  At  this  station  the  two  zone  model  aits  as  i 


one-zone  model. 


It  is  possible,  however,  to  introduce  a  non mst  ant 


32 


value  of  the  discriminating  scalar  such  that  both  zones  always  contain  a 
non-negl ibi le  amount  of  energy  as  shown  in  Appendix^.  Figures  14a-c  show 
the  development  of  the  probability  forg<S,  calculated  as  the  integral 
from  £*C>  to  €■*£,  of  the  log-normal  distribution.  Due  to  the  decrease 
in^Othis  probability  increases  downstream.  In  fig.  15  the  interzonal 
energy  transfer  is  shown.  Tnis  energy  transfer  is  negative  over  the  wnole 
cross-section  (x/0  »  60)  which  means  that  energy  is  transferred  from  zone 
1  to  zone  2,  but  it  is  not  strong  enough  to  keep  up  with  decrease  of 
due  to  the  decrease  in  the  dissipation  rate.  Figure  16  shows  the  mean  and 
the  standard  deviation  of  the  dissipation  rate.  The  stanoard  deviation  is 
everywhere  larger  than  the  mean  value  as  calculated  from  the  refined 
similarity  hypothesis  (eq.  (4.17)). 

Cone lusions 

A  new  multi-scale  (multi-zone)  concept  for  turbulent  flows  nas  been 
developed,  and  exact  equations  valid  inside  the  different  zones  are 
derived  ana  some  simple  initial  modelling  assumptions  have  been 

constructed  in  order  to  close  the  system  of  equations.  The  model  i  ■> 
applied  to  calculate  decaying  isotropic  turbulence  and  a  plant  jet.  T . 
results  are  qualitatively  in  aqreanert  with  expectants,  but  due  ti  the 

complete  lac*  of  experimental  data  it  is  not  possible  t  Ira*  an,  further 

conclusions.  More  wo  r*  is  probably  neiessary  i n  order  tr  refine  tne 

mode  I  i  l  ng  and  to  e  *  t  png  the  mu  1 1  i  -  s.  a  t  e  i  ont  ep  t  t  <  >  ,  ec  u  ■>;  roe  f  »-*  1 

Ac  know  1 edgenent 

Ini',  *iir«  was  supported  Uy  •cluse  .,r  vit  \  i . 
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APPENDIX  1 


PDF  equation  for  c 

Let  £  be  given  as  .  An  exact  transport  equation 

for  can  be  derived  Dy  taking  the  derivative  of  the  instantaneous 
momentum  equation 

w.r.t.  Mj.  and  multiply  the  result  with  te  ■ 

The  result  is 

4?  *  V.4?  -  -  2^4,14  g/° 

~2v&i  n) 

where  the  pressure  can  be  related  to  the  velocity  field  by  the  Poisson 
i ntegra 1 . 


Oenotinq  by  £  and  ^  the  probability  space  variables  corresponding  to  £ 
and  ,  respectively  and  defining  a  fine-grained  joint  Pdf  for  velocity 
and  dissipation  as 

3 

P £,*■,*)  s  tXf?'  -  F )  //  cfy&L  -  ^ 

the  equation  tor  *,  i  is  easily  obtained 

by  jsmg  standard  operations  with  tne  #^-f  unt  1 1  on  . 


-S -S -V  -t  -S  -V  -  -  --V -  ■  Al  .l~.iv-V  -V-  ■  -  ■ 


3' 


3tp  +  sin^py 


- »£ ae&py  - ^-£ 

*C 

-y&fjfixvi  &&$]  +2v&  rffiafyo 

% 

+  v£t&%.%  *  vfahV'Zv,)?’ )  (2) 


By  integrating  this  equation  over  the  velocity  space  the  equation  for  the 
marginal  Pdf  of  £  ,  is  obtained  as 

•2**  *  W*S£>P)  "  ~  *$ 

+  2vJe  %/>&*.  *  +  v %vAv  &,) p y 

where 

5  s  Jyttij)  -  e  y 


APPENDIX  2 


Analytical  solution  for  the  variance  of  g.  in  the  case  of  isotropic 
turbulence 

Consider  the  following  form  of  the  equation  for  the  variance  of  &  in 
the  case  of  isotropic  turbulence. 


df£j> _ c,  . .  t>  fl) 

di  t  y 

where  2^^^-^and  where  and  ('S)  can  be  found  from  thef*)  -  <£,> 

model  as 


f*> 


a) 


<e> 

(e>c 


‘  ft  + 


Cc 


wnere  ~  ^ O) 


and 


b 


•  SOo 

The  solution  of  (1)  is  then 


Gn>o 


{*) 


From  (2)  and  (3)  it  can  be  seen  that  the  turbulent  Reynolds  number 
r><2%>dj>decreases  with  time  when  1.0  <"  <  2.0.  The  commonly  used 

value  is  *  1.9  which  is  consistent  with  a  decrease  of  Re.  The 
requirement  that  the  ratio  )  is  increasing  with  increasing  Ke 

gives  the  following  condition  on  C. 


Cj  >  2CZ 


f-5  ) 


leading  to  the  following  condition  on  A  from  (4.11) 

A  >  O 


The  constant  can  also  be  related  to  the  constant yx  in  the  Kolmogorov 
expression  for  the  variance  of  £  by  the  requirement  that  the  ratio 
>  calculated  from  (3)  and  (4)  should  have  the  same  time 
dependence  as  (4.17).  The  result  is 


c$  -  2ci  (*+  TS*  *0*  ) 
The  constant  A  in  (4.11)  can  then  be  given  by 


(7) 


A  -  -yyU  -  Cz  )  fg) 

The  constant y U  can  further  be  related  to  the  fractal  dimension  of 
intermittent  turbulence. 


APPENDIX  3 

Integrals  and  moments  of  the  loq-normal  distribution 


Let  £.  be  a  stochastic,  log-riormally  distributed  variable, 
logarithm  of  £,  is  then  normal  distributed  and  the  Pdf  of 
is  given  by 


'  -rH^rs 


Ae  TFir 


ATT 


where 


m*  ,  Az  “  { A&+.S  —  £ >J*> 


The  Pdf  of  E y  PtK)  can  then  be  determined  os 

p*f‘>  -  *  p<*> 


which  gives  for  P(£) 


P(t>  - 


A  eyTr 


Consider  the  integral  X' 

€./ 


T6'  2  J  a/e  ep  Pre) 


Introducing  (2)  gives 


TA  '  aAfA  r^e -"A 


Introducing  a  new  variable  U(£)  defined  by 


Ufe) 


lece  -  nq 

A7T 


c/u 


c/£ 

AeTT 


gives  the  following  expression  -for 


ufe,) 


~y~=r- ■  J  o/u  APt->  ~  °2‘) 


—  04 


This  integral  can  be  further  transformed  into  (see  Ref.  [15],  p.  303) 


e*y>  +  £/3>z/>z)j[e+J£{o(et)  —  ^  )  +/] 


where  )is  the  error  function  defined  as 

v 

*  wfdi  e~iZ 

o 

The  moments  of  £  can  now  be  calculated  as 


<ep>  -  T'’ 

60 


+  jpzAz 

e 


/*) 


By  combining  the  first  and  second  order  moments  the  parameters  m  an 
can  be  expressed  as 


A 


Z 


U  f  J 


<€*>  x 

<e> *■  ) 


rn.  *  -  (p 

7/7zm> 

/  ay 


n 


where  /£>*)  „  _  &))*>  . 


The  integral  _ZT^  must  be  calculated  numerically,  this  can  easily  be 
done  by  one  of  the  approximations  for  the  error-function  given  in  Ref. 
15.  In  this  work  the  following  approximation  is  used. 

z  5- 

-  /  -  e  *  £cs  l*  y  >o 

where  if*  J+px  »  an<^  constants  are 

P  =  0.3275911  ,  ax  =  0.254829592  ,  a2  =  -  0.284496736 
*3  =  1.421413741  ,  a4  =  -1.453152027  ,  a$  =  1.061405429 
For  y<0  the  symmetry  properties  of  can  be  used 

&+jL(y)=  —  e+fi 


APPENDIX  4 


Transformation  of  the  multiscale  equation 

By  defining  A  *  kM,Pfe)  the  modelled  equation  for  zone  conditioned 
kinetic  energy  can  be  written  as 


^4  y-  ^>^4  -  ^ 


n) 


where  is  the  collection  of  production/dissipation  and  interzona' 

transport  terms. 

A 

A  new  discriminating  scalar  £  defined  as 


£  -  £ 


* 

and  a  new  dependent  variable  4  defined  as 


A  «  k'"Pf,i)  ,  =  <e>Pg(i) 


are  introduced  such  that 


JoteAfe)  *  {/<) 


> i  ^ 

Then  fol lows  4  *  and  the  equation  for  A  becomes 


a) 


where  ^r>  is  the  source  term  in  the  -equation . 


t-juit-un  U  ;  '  •>  now  tr4nsf  ..nnej  tc  a  ne*  oor^mate  svster  :  y  r 


t rms f  ormat  ’  jn 


^  ,*.e  1 


°  ,  e 


-r*re  K  m  **  m  *  ^  ~  <k> 


js t ne  t;  ■  )»  ■  og  •  racist  '-T.at  1  in  •■  e>, 


0  f  1  ° 

“"«  *  <•  o7>  %  ^ 


•i  '4  ~  *5  ^>4 


t.ne  f  j  i  ijn'i^j  trin^t  jnieo  «*\H»at'.,  n  '•ej  v . » 

^  1  ~  2 '%>£<*>  if 


'7f>‘  Z  *  '/'a  )  +&)$r>  ? //tJ!  1  +  Z 


•-  *  *r‘-  ’->»  >*■"*•■■:.  n*.  ♦  r  f  <  )  ar  :  ,,(>  *  '  E  >  I 


’’  ^  *'  no»  tne  *.juat  •  '  ♦.)«•  ♦  n*.  «  •  net 


►"v^y  ’  f  ’one 


.  :  •  a  -  h  ;  "1^  -  j„  a*  ■  ,»  *  ,.-  k 


•  r-  , .  , 


'■'!<?„»■  !t  ■  nj  'yf*-  *  ir  •  n  t  e  r  « a  ^ 


^  >4  *■  *  <*  ^  -  2  &*>  4  <*>*«£*]„ 


+  7^1  <*>*('&'>  ~:a  +  r  L  A  J 


If  an  infinite  number  of  zones  is  used  the  form  of  Af£)  is  known,  but  in 

A 

the  cas  of  a  finite  number  of  zones  A(()  must  be  approximated. 

when  only  two  zones  are  used  we  suggest  the  following  approximation 
f  or  )  . 


A 

A  “ 


2  kj  F  for  OFF*/ 

lilO-ir)  for/*£<EK 


A' 


*  £  .  ,  kt 

where  is  given  by  m  /  ~r 

Tms  approximation  is  illustrated  in  fig.  1, 


wmcn  gives  »  2  k, 

/«/ 


c 

*  f  Je  *  / 


an  be  approximated  as 


~  z  / y  J "  k  f ) 
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Fig.  12  Evolution  of  P (€,)  as  calculated  from  the  transport  equation 
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